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Chapter 1 
Introduction 



A wide range of problems in modern high energy physics and cosmology involves 
the dynamics of quantum fields in highly out-of-equilibrium conditions, includ- 
ing relativistic heavy ion collisions [1-3], quantum fluctuations in inflationary 
cosmology [4-6], preheating after inflation [6], and baryogenesis [7]. It turns 
out that the traditional methods of (vacuum or thermal) quantum fleld theory 
(QFT) are not suited to describe these complicated non-thermal processes, and 
a new framework of nonequilibrium QFT [8,9] is needed. Even though a knowl- 
edge of the complete dynamics of interacting quantum fields is clearly beyond 
reach, there is a variety of approximative methods that are able catch the es- 
sentials of the problems under study. In a so called kinetic regime with weak 
interactions and slowly varying classical backgrounds the standard methods of 
quantum kinetic theory reduce the problem considerably to solving the famous 
(quantum) Boltzmann transport equations. For many problems of interest these 
transport equations provide a remarkably good approximation for the essentials 
of nonequilibrium quantum dynamics. However, certain problems are inherently 
very sensitive to quantum coherence (or interference), quantum reflection from a 
potential being a typical example. The Boltzmann equation approach inevitably 
loses the effects of nonlocal quantum coherence, and thus is not very well suited 
to study for example the quantum reflection problem in electroweak baryogenesis 
or the particle production in preheating. 

In this thesis we present a novel approximation scheme related to the quan- 
tum kinetic theory, that enables us to treat nonlocal quantum coherence in the 
presence of decohering collisions with simple enough Boltzmannian-type trans- 
port equations. The key element in our scheme is the flnding of new singular 
shell solutions in the phase space of 2-point correlation function, that are lo- 
cated at fco = foi' spatially homogeneous problems and at = for a static 
planar symmetric case. When the complete phase space structure, including 
these new coherence solutions in addition to the standard mass-shell contribu- 
tion, is inserted in the Kadanoff-Baym (KB) equations for the correlator, we 
obtain a closed set of transport equations for the corresponding on-shell distri- 
bution functions, thus giving an extension to the standard quantum Boltzmann 
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equation to include nonlocal coherence. 

The thesis consists of three original research papers [I]- [III] and an intro- 
ductory and summary part presented below. In chapter [2] we introduce the 
basic mechanisms of (electroweak) baryogenesis and preheating, that are good 
examples of highly nonequilibrium processes in the early universe. In chap- 
ter [H] we present the basic formalism of nonequilibrium QFT needed to derive 
the KB-equations for two point correlation functions, using the two-particle ir- 
reducible (2PI) effective action method. Chapter H] then presents a detailed 
survey of the main contents of our work, the novel approximation scheme. In 
chapter El we review a few applications that we have so far considered using 
our formalism, including the Klein problem, (coUisionless) quantum reflection 
from a CP- violating mass wall, and examples of coherent production of decaying 
fermionic and scalar particles relevant for preheating. Finally, chapter [6] contains 
conclusions and outlook. 
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Chapter 2 

Nonequilibrium processes in the 
early universe 

According to modern theories of cosmology and particle physics the expand- 
ing universe has once been in an extremely dense and hot state (Hot Big Bang 
scenario), consisting of quantum plasma, which during the major part of the 
early evolution is very close to thermal equilibrium. Besides this overall picture 
of thermal plasma, however, many crucial processes in the early universe are 
inherently highly non-thermal, including inflation, preheating, and baryogene- 
sis. The careful understanding of these processes is of primary importance in 
modern cosmology, providing an important field of applications to the methods 
of nonequilibrium quantum field theory. In this chapter we introduce the basic 
mechanisms of baryogenesis, focusing on a model called electroweak baryogenesis 
(EWBG), and the process of preheating after inflation. 

2.1 Baryogenesis 

The visible matter content of the universe, such as planets, stars and interstellar 
gas, consists of protons, neutrons and electrons. In astrophysics it is classified 
as baryonic matter, since the bulk of the mass is in protons and neutrons that 
are baryons. There is strong evidence that no large domains of antimatter exists 
in the universe [10, 11], implying that the universe has an excess of baryons 
compared to antibaryons which is called baryon asymmetry. The combination 
of data including several experiments of the fluctuations of cosmic microwave 
background (CMB) gives the following experimental measure of this asymmetry, 
the average baryon to photon number ratio in the universe [12]: 

^ = (6.ltg:^) xlO-i°. (2.1) 

Of course it might be possible that the baryon asymmetry is an initial condition 
in the evolution of the universe. Despite being very unnatural, this explanation 
for the baryon asymmetry is not consistent with the cosmological inflation [4,5], 
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which is one of the backbones of modern cosmology explaining the homogeneity 
and flatness of the universe as well as the primordial density fluctuations that 
will give rise to structure formation and the observed fluctuations in the CMB 
spectrum. The problem is that the exponential increase in the size of the universe 
by at least a factor e^'^ during the inflationary period dilutes any prior baryon 
number to totally negligible level. For this reason, we are very tempted to seek 
out different ways of creating the baryon asymmetry in the universe after the 
period of inflation. 

A process that gives rise to a permanent baryon asymmetry at cosmological 
scales is called baryogenesis. The idea of such a process originates from Sakharov 
[13], who presented three conditions that any model for baryogenesis should 
necessarily fulflU: 

1. Baryon number violation. 

2. C and CP symmetry violations. 

3. Departure from thermal equilibrium. 

The first condition is obvious. If the second is not fulfilled, then for every reaction 
producing particles there is a counter-reaction that produces antiparticlcs at the 
same rate. The third condition is the most interesting for the scope of this work. 
It follows from the CPT-theorem that the masses of particles and antiparticles 
are equal, and consequently the thermal average of the baryon number will 
vanish in equilibrium. We conclude that every scenario for baryogenesis must 
be a nonequilibrium process. 

Several models for baryogenesis have been proposed that fulfill the Sakharov 
conditions (for a recent review see e.g. . [7]), including GUT baryogenesis [4], 
electroweak baryogenesis [14], leptogenesis [15] and Afflect-Dine baryogenesis 
[16] as the most prominent candidates. The first of these, GUT baryogenesis, is 
based on decays of heavy gauge bosons with masses of order Mqut ~ 10^^ GeV. 
While it provides a scenario fulfilling all the Sakharov conditions, it has serious 
problems with inflationary models related to the high reheating temperature 
required, and the consequent overproduction of gravitinos [7]. The latter mod- 
els, electroweak baryogenesis, leptogenesis and some variants of Affleck-Dine 
baryogenesis are based (directly or indirectly) on electroweak baryon number 
violation [17], which is a quantum anomaly in the electroweak sector of the 
standard model allowing the baryon number to be badly violated at high tem- 
peratures. These models differ however substantially in the mechanisms of how 
the required out-of-equilibrium conditions are reached. In what follows we will 
focus on EWBG in more detail, trying to elaborate the basic mechanism and 
the necessity to use the methods of nonequilibrium quantum field theory in its 
study. 
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2.1.1 Electroweak baryogenesis 
Electroweak baryon number violation 

Baryon and lepton numbers are classically conserved in the standard model. 
However, at the quantum level this is not the case. It can be shown that the 
axial current in the electroweak sector of the standard model and consequently 
the total baryon and lepton number currents and are anomalous i.e. not 
exactly conserved [18,19]: 

d,fB = d,fL = Nf [^.K^W^'^'' - ^2^^-^^"') ' (2-2) 

where Nj is the number of fermionic families, PV^^, and F^,^ are the field strength 
tensors of the SU{2)l and f/(l)y gauge symmetries with the duals W""'^" = 
Ig/ii/po-pyo^ and an analogous expression for F, and g and g' are the associated 
coupling constants, respectively. Equation (12.21) implies that the total change in 
baryon (lepton) number from time t = to some arbitrary final time t/ is given 
by: 

AB = AL = Nf [Ncsitf) - iVcs(O)] - Nf [ricsitf) - ncsiO)] , (2.3) 

where 

^cs = d'xe^^^ {w^^Al - ^ge,,,A1A)Al^ 

^cs = d'xe'^'F,,Bk (2.4) 

are called the Chern-Simons numbers of the SU{2)i and f/(l)y gauge symme- 
tries. We proceed by considering the vacuum structure of the gauge fields 
and B^. It turns out that the abelian U{1)y sector has a trivial nondegener- 
ate vacuum with 5 = 0, but the SU{2)l sector instead has a discrete set of 
degenerate vacua with 

A=-g-^Vgn, (2.5) 

where gn{x) = e^^'fi^)^-'^/'^^ n is an integer and Tj are the generators of the SU{2) 
gauge group. Using the vacuum structure Eq. (12.51) in Eqs. (I2.3I) - (I2.4I) we see 
that the change in baryon number in transitions between the different vacua is 
given by 

AB = NfANcs = Nf^Jd'xe'^''TT[g-%g^g;;'d,g,,g-'dkgn\ 

= NfU. (2.6) 

We see that the baryon number changes in integer multiples of Nf in the 
transitions between different vacua. The structure of the effective potential 
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Figure 2.1: The vacuum structure of the SU{2)l gauge fields. The set of de- 
generate vacua can be labelled with an integer Chern-Simons number Ncs- At 
zero temperature the transitions between different vacua by quantum tunneling 
through the potential barrier are extremely suppressed. At high temperatures 
the transitions are possible by thermal activation via sphaleron configurations. 

for the gauge fields is sketched in Fig. 12. 1^ where the minima correspond to 
different vacuum configurations labelled by the (integer) Chern-Simons number 
Ncs- But how could the transitions between different vacua actually take place? 
At zero temperature the only possibility is by quantum tunneling through the 
potential barrier. This corresponds to the so called instanton configuration [17], 
but it turns out that the tunneling rate is neglig 

With this rate not a single proton could have been produced in the lifetime of the 
universe! At high temperature the situation is better. The transitions between 
the vacua can take place through a thermal activation over the potential barriers, 
via the so called sphaleron field configuration [20,21]. The thermal transition 
rate corresponding to this process is shown to be [14, 22] 

Fsp ~ T^e-^-(^)/^, (2.7) 

where E^p is the energy of the sphaleron configuration, which is related to the 
vacuum expectation value of the Higgs field (0) by E^p/T ^ 4O(0)/T. This result 
is valid only in the broken phase i.e. when the gauge symmetry is spontaneously 
broken and the gauge bosons are massive. In the symmetric phase with massless 
gauge bosons the sphaleron rate is instead given by (in the minimal standard 
model) [23] 

Fsp^ (25.4±2.0XT^ (2.8) 
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where ~ 1/30 is the weak couphng constant. To see if the sphaleron transi- 
tions are fast at the time scales of the expanding early universe, the sphaleron 
rate of a unit comoving volume: Fsp/T^, needs to be compared with the Hubble 
expansion rate of the (radiation dominated) universe [4] 

H = lMgl/'—, (2.9) 
mpi 

where 5'*(T) counts the total number of effectively massless degrees of freedom 
and mpi = 1.22 x 10^^ GeV is the Planck mass. In the minimal standard model 
at high temperatures of order T > 100 GeV we have g^: = 106.75 so that the 
symmetric phase sphaleron rate in Eq. (12.81) is very large compared to the Hubble 
rate. Moreover, if {4>)/T < 0{1) then also the broken phase sphaleron rate in 
Eq. (12. 7p is large, and the first Sakharov condition is fulfilled in both of the 
phases. Actually, it is crucial in EWBG that the broken phase sphaleron rate 
is smaller than the Hubble rate so that the the generated baryon asymmetry is 
not washed out. Next we will briefly consider the electroweak phase transition 
between the symmetric and broken phases, where the gauge bosons and fermions 
become dynamically massive. This transition takes place at the (electroweak) 
temperature scale T ~ 100 GeV, and it provides a scheme for the creation of a 
permanent baryon asymmetry through EWBG, if the transition is of first order. 

Electroweak phase transition 

In the electroweak theory of the standard model [24-26] the masses of the gauge 
bosons and Z° and all fermions are generated through spontaneous symme- 
try breaking of the gauge symmetry SU{2)l x f/(l)y. The "classical" potential 
of the {SU (2) L doublet) Higgs scalar field 

n0) = -^0^ + ^0\ (2.10) 

where (f) = v2$t$ ^ is minimized for (j) = vq = fi/ \/X corresponding to a degen- 
erate set of vacuum configurations. Choosing a particular vacuum from this set 
breaks the gauge symmetry spontaneously giving rise to masses for the gauge 
bosons, which are proportional to the vacuum expectation value vq of the Higgs 
field. This is how the Higgs mechanism works at zero temperature. 

At finite temperatures the "classical" Higgs potential Eq. (I2.10p gets tem- 
perature dependent quantum corrections that will become more and more im- 
portant when the temperature increases. These corrections are taken in account 
by computing the free energy (or the effective potential) of the Higgs field [7]: 

VM, T) = D{T^ - T^)<i>'' - ET<j)' + + ... (2.11) 

where the values of the parameters To, D and E depend on the considered model. 
We have only given the dominant terms from perturbative calculations. The 
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Figure 2.2: The Higgs field effective potential for a first order phase transition 
at different temperatures. Because of the barrier between the local minima, the 
phase transition takes place at a temperature T < releasing latent heat. 

temperature dependence of the potential Eq. (12. lip gives rise to an interesting 
behaviour of the Higgs field. Let us first consider the case with vanishing third 
order term E = 0: when T > Tq the second order term is positive implying that 
the minimum of the potential and the corresponding vacuum expectation value 
(VEV) of the Higgs field is zero: (0) = 0, while for T < Tq it is finite: (0) ^ 0, 
since the negative second order term dominates for small 0. At the critical tem- 
perature T = Tc = Tq there will be a phase transition from the former symmetric 
phase to the latter broken phase. This is called electroweak phase transition, and 
for the case ii^ = it is a second order transition with smoothly increasing VEV 
(0) as the temperature decreases. The case with E > presented in Fig. 12.21 
is more interesting for us. At the critical temperature T = Tc there is now a 
potential barrier between the symmetric and broken minima. Because of this 
barrier the phase transition actually takes place at a lower temperature T < T^, 
and latent heat is released due to the energy difference between the minima. 
Hence the transition is of first order and it happens by nucleation of broken 
phase bubbles inside the symmetric phase bulk. These bubbles start to grow 
rapidly reaching soon a stationary expansion speed. This scenario with first 
order phase transition is of crucial importance for the electroweak baryogenesis. 
It is just around the edge (or wall) of these broken phase bubbles where the last 
two Sakharov conditions of CP violation and departure from thermal equilib- 
rium are fulfilled. Unfortunately, it has been shown by nonperturbative lattice 
calculations that for the minimal standard model the phase transition is not of 
first order [27]. For this and other reasons [e.g. not enough CP violation) we 
have to seek new possibilities for electroweak baryogenesis in extensions of the 
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Figure 2.3: Sphaleron rate (red), vacuum expectation value of the Higgs field 
(green), and the source asymmetry (yellow) in the bubble wall region as a func- 
tion of the coordinate normal to the wall front. 



standard model, such as the minimal supersymmetric standard model (MSSM). 
Next we will give a brief conceptual description of the mechanism of (eletroweak) 
baryogenesis for a generic model with a strong first order phase transition. 



Description of the EWBG mechanism 

In Fig. I2.3l we present a schematic cross-section from the expanding bubble wall 
front where baryogenesis takes place. In the bubble wall region the VEV of the 
Higgs field and consequently the masses of the fermions (green line) are spatially 
dependent changing smoothly from zero in the symmetric phase to a finite value 
in the broken phase. The varying complex mass of a fermionic eigenstate gives 
rise to the required CP violating effects that will generate an asymmetry in 
left chiral number densities between particles and antiparticles (yellow line). 
This source asymmetry then creates a pseudo chemical potential that biases the 
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sphaleron transitions in front of the wall in the symmetric phase to produce more 
baryons than antibaryons. This net baryon number is not washed away in the 
broken phase, because the sphaleron rate (red line) is too small there compared 
to the Hubble expansion rate H. At the end, when the bubble wall has passed 
and the plasma is back in thermal equilibrium, a net baryon number density has 
been created. 

One of the most difficult problems in the actual calculations of the baryon 
asymmetry is to find out the source asymmetry due to the CP violating ef- 
fects. An accurate calculation would necessarily involve the use of nonequilib- 
rium quantum field theory, and this is the background motivation for our work. 
In earlier works the problem has been studied in the semiclassical WKB ap- 
proach [28-33] and later with the methods of quantum kinetic theory in [34-38], 
both approaches using the Boltzmann transport equations for CP- violating 
phase space densities. These methods should provide a solid approximation 
in the semiclassical limit i.e. the case of a thick wall compared to the mean 
free path of the interacting fermions. However, in the thin wall limit the dom- 
inant source for the asymmetry comes from the quantum refiection processes, 
which are inherently nonlocal and absent in the standard WKB and kinetic 
approaches. Attempts have been made to treat the refiection phenomena by 
including collisions in the Dirac equation [39-43], but no consistent framework 
based on quantum field theory has been introduced. In this work we present an 
approximation scheme based on the kinetic approach that enables us to treat the 
quantum refiection in a simple but consistent way in the presence of decohering 
collisions. 



2.2 Preheating 

Cosmological infiation [4,5,44-52] (for a recent review see e.g. [6]) is a period 
of rapid exponential expansion in the very early universe during which the size 
of the universe increases by a huge factor (at least e^°). As mentioned in the 
beginning of last section, it provides a natural explanation for the homogeneity 
(horizon problem) and fiatness of the universe. It also explains the (almost) 
scale invariant primordial density fiuctuations that will give rise to large-scale 
structure formation and the observed anisotropies in the CMB spectrum [12]. 
Most of the infiationary models are based on the peculiar dynamics of one or 
several scalar fields, called infiaton(s), whose vacuum condensate dominates the 
evolution of the universe in a so called "slow roll" phase, causing the exponential 
growth. Because of the huge and very rapid expansion the universe is typicall}0 
in a highly non-thermal and very cold state at the end of infiation. However, 
the baryogenesis scenarios require energies greater than the electroweak scale 
and the primordial nucleosynthesis requires that the universe is close to thermal 

^The warm inflation scenario with particle production during inflation is an exception 
[53,54]. 
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equilibrium at the temperature ~ 1 MeV at some stage after inflation. A mech- 
anism for reheating the universe after inflation is thus needed to retain the Hot 
Big Bang scenario. 

The modern scenarios of reheating consist of a preheating stage followed by 
thermalization. In the preheating stage [55-60] the inflaton condensate goes 
through rapid oscillations, so that the couplings to other matter fields give rise 
to particle production via parametric resonance. The basic mechanism is simple: 
the coupling to the infiaton(s) gives rise to rapidly oscillating effective masses 
for the matter fields that will bump up the particle numbers exponentially for 
certain momentum modes in close analogy with the Floquet theory of growing 
exponents [61,62]. The modes in these resonance bands will quickly obtain huge 
occupation numbers (for scalars with no Fermi blocking). This rapid growth of 
perturbations is followed by backreaction and rescatterings. Backreaction means 
the effects of the growing perturbations (particle production) back on the dy- 
namics of the inflaton condensate. In most scenarios it will rather quickly shut off 
the inflaton oscillations (faster than the Hubble expansion rate of the universe) 
and consequently terminate the particle production. Before that the fluctuations 
of the inflaton field itself will grow and give rise to rescatterings i.e. couplings 
between different momentum modes leading to the growth of occupation num- 
bers for non-resonance modes as well. When the oscillations of inflaton field shut 
off completely the preheating stage ends and the fields start to thermalize. This 
thermalization process can be very complex, including regimes of driven and free 
turbulence [63], which makes it difficult to estimate the final reheat temperature. 
Moreover, in certain cases the universe might enter a "quasi-thermal" phase with 
a kinetic equilibrium reached much before the full chemical equilibrium [64]. 

In addition to scalar particles, the parametric resonance during preheating 
may produce a significant amount of fermionic particles. This resonant produc- 
tion could possibly lead to dangerous relic abundances of problematic particles 
such as gravitinos [65,66]. Fortunately, extensive studies have shown that grav- 
itino over-production can be avoided during the preheating in realistic supersym- 
metric theories [67-73]. Another interesting aspect in the fermionic preheating 
is the possibility to generate heavy fermions with masses of order 10^^-10^^ GeV, 
that could be important in e.g. leptogenesis [74,75]. In chapter [5] we will apply 
our approximation scheme with decohering interactions to study a simple model 
of fermionic preheating, during which the fermion is subjected to decays. 



11 



12 



Chapter 3 



Basic formalism of quantum 
transport theory 

The standard methods of (vacuum or thermal) quantum field theory (QFT) are 
not well suited for the study of nonequilibrium quantum fields for several rea- 
sons. First of all, the basic quantities of interest in nonequilibrium QFT are 
expectation values of operators in contrast to transition amplitudes in the stan- 
dard vacuum QFT, requiring the formal extension of the time variable into a 
closed time path with two different branches. Other specific issues in nonequi- 
librium QFT are related to for example secularitj0 causing the complete failure 
of the standard perturbative expansion in many cases of interest [76-78] . In this 
chapter we introduce the basic concepts of nonequilibrium quantum field theory 
in order to derive the Kadanoff-Baym (KB) transport equations for fermions 
and scalar bosons. We start by introducing the closed time path formalism with 
four basic propagators. Then, we use the two-particle irreducible (2PI) effective 
action methods to derive the self-consistent (Schwinger-Dyson) equations of mo- 
tion for the full 2-point correlation functions of the system, which we then write 
in the form of KB-equations. Finally, we list some physical observables that can 
be expressed in terms of these 2-point functions. 

3.1 Closed time path formalism 

The closed time path (CTP or Schwinger-Keldysh) formalism was developed by 
Schwinger [79] and Keldysh [80] and further refined by many infiuential works, 
including [81-89]. The basic idea of the formalism is simple: In order to study the 
expectation values instead of transition amplitudes by the methods of quantum 
field theory, the time coordinate must be extended to a closed time path from 

^The standard perturbative approach suffers from secular terras giving rise to big late- 
time contributions in all orders of perturbation expansion, no matter how small the coupling 
constant is [76-78] . Heuristically, this can be seen in the sense that et ^ 1 for a big enough t 
for every e (no matter how small). 
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Figure 3.1: The closed time path of the CTP formalism. 

initial time to to final time tj (often taken to be oo) and then back to (see 
Fig. 13. ip . The need for this closed time path can be demonstrated by writing 
the expectation value of a real scalar field (f){x) in a state defined by arbitrary 
density operator p in terms of a path integral (we choose here = 0): 

(0(x)) = Tr{0(x)p} = Tr{f/(t;,O)?7(O,t)0s(x)?7(t,O)ps(O)f/(O,ty)} 



#/rf0jd0o#'o[(0/|t/(t/,t)|0t)0(x)((/)t|[/(t,O)|0o) 
x{<l>o\psm<l>o){<t>o\UiO,tf)\<t>f) 



d(j)f d(f)o d(f)Q p[0o(5),0o(x)] 



4>o{x) 
4>f{x) 



0^(x) expliSf^^]} 



X / D(j)- ex^{iS[(p-]*] 

l^'^{x) 

D0+D0-p[0+(O,f),0-(O,f)] 0+(x)exp{i(^[0+] 

(3.1) 

where p[0^(O, a;), 0~(O, x)] = {(j)'^\ps{Q)\4>~) is the initial density matrix and 
the subscript Os denotes an operator in the Schrodinger picture in contrast 
to the Heisenberg picture without a subscript. We see that the path integral 
representation involves two "histories" , for which the evolution is chronological 
from to and antichronological from to 0, respectively. The field values in 
these +/— branches are independent except the boundary condition 0+(t/,x) = 
4>~(tf, x) closing the path at t = tf, which we did not write explicitly in the last 
row of Eq. (l3All . 



^For this reason the CTP formahsm is often called "in-in" formalism on the contrary to 
the traditional "in-out" formalism of quantum field theory with transition amplitudes between 
incoming and outgoing states. 
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3.1.1 Propagators 

It turns out that the higher n-point Green's functions of the system will auto- 
matically become time ordered along the closed time path. Of special interest 
in quantum field theory are the 2-point functions or propagators, which are now 
defined as 

i^cM = (Tc ^TT{pTc[<Piu)<Piv)]} (3.2) 

iGc{u,v) = {Tc[^{u)i^{v)]) = Ti{pTc[^{u)i,{v)]} (3.3) 

for a real scalar field and a fermionic field ijj, respectively, p is some unknown 
quantum density operator describing the statistical properties of the system, 
and Tc defines the time ordering along the closed time path C, shown in Fig. 
13. H in the sense that the points on the lower (negative) branch are "later" than 
those on the upper (positive) branch. When written in terms of the ordinary 
real time variable running from — oo to +oo, the closed time path propagators 
f l3.2p -( !373l) contain four distinct contributions depending on the time branches 
of the "complex" CTP-coordinates u and v. That is, using indices a, 6 = ± to 
label the positive/negative branches, the scalar propagators are decomposed as: 

iA+"(M,f) = iA<(M,f) = {(f){v)(f){u)) 
iA~^{u,v) =iA>{u,v) = {(f){u)4){v)) 

iA'^^{u,v) = iAF{u,v) = 9{uo — vo)iA^{u,v) -\- 9{vq — uo)iA'^{u,v) 
iA {u,v) = lApiu^v) = 9{vo — uo)iA^{u,v) + 9{uo — vo)iA^{u,v) , 

(3.4) 

where now and are ordinary time coordinates. Similarly, for fermionic 
propagators we have: 

iG+~{u,v) = -iG<{u,v) = -{iIj{v)iIj{u)) 
iG-+{u,v)= tG>{u,v) = {iIj{u)^{v)) 

iG^^{u,v) = iGp{u,v) = 9{uq — vq)G^{u,v) — 9{vq — uq)G'^{u,v) 
iG {u,v) = iGp{u,v) = 9{vo — uo)G^{u,v) — 9{uo — vo)G'^{u,v) . 

(3.5) 

Using the generic notation Q = {G,A} to denote fermionic/scalar popagators 
we see that Qp and Qp are the chronological (Feynman) and anti-chronological 
(anti-Feynman) propagators, respectively, while and are the so called 
Wightman functions. In our further analysis we are especially interested in the 
dynamics of these Wightman functions, which contain the essential thermal or 
out-of-equilibrium statistical information of the quantum system under study, 
in order to compute for example the expectation values of the number current 
and the energy momentum tensor T^'^ . 

Before we start to build the calculational scheme of the OTP formalism, let 
us introduce a few more Green's functions, which are useful in the following 
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analysis, and list some of their properties. First, we define the retarted and 
advanced propagators: 

Q\u,v) = Qf-G^ = -e{v''~u^){g>TG^), (3.6) 

where now =]= refers to bosons/fermions. The definitions 03.41) . fl3.5p and 03.61) 
then imply that the propagators have the following hermiticity properties: 

[tA<'>{u,v)]^ = iA<'>{v,u) 
[^G<'>(M,t;)7°]^ = «G<'>(t;,M)7°, (3.7) 

and further [iA'^{u,v)f = —iA"'{v,u) and [iG^{u,v)'y^]^ = — zC(f,u)7°. These 
latter identities for retarted and advanced propagators suggests to decompose 
them into hermitian and antihermitian parts: 

Gh = liG^' + Gn 

A = ^^{G''-G1 = 1{G>tG<). (3.8) 

The antihermitian part A is called the spectral function. Based on Eqs. 03. 6p it 
is easy to show that Gh and A obey the spectral relation: 

Gh{u, v) = -isgn(M° - v^)A{u, v) . (3.9) 

3.2 Two-particle irreducible effective action 
and Schwinger-Dyson equations 

In a nonlinear quantum field theory, including the majority of the interact- 
ing theories, the 2-point correlation functions necessarily couple to higher order 
correlators and so on, to form an infinite Schwinger-Dyson hierarchy of equa- 
tions analogous to Bogoliubov-Born-Green-Kirkwood-Yvon (BBGKY) hierarchy 
in classical statistical mechanics [9]. For this reason it is an impossible task to 
solve the 2-point correlators exactly; that would correspond to a full solution 
of the nonlinear QFT. A major paradigm in practical applications thus is to 
truncate this hierarchy (slave the higher order correlators) in an appropriate 
way. The truncation can be done in many ways, for example by a brute use 
of perturbation theory or a loop expansion. However, these standard methods 
do not generally provide a good approximation for out-of-equilibrium dynamics, 
because of several problems e.g. with secularity [78]. 

A way to evade these problems is a method of obtaining the (truncated) 
equations of motion from variational principles of increasing complexity. On the 
first level we obtain an equation of motion for the field expectation value (0) 
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only, in the next level to (0) and 2-point function Q ~ {{<P — (0))^)? and so 
forth. The effective action corresponding to the n-th level of this hierarchy is 
called n-particle irreducible (nPI) effective action Fnpi. The heuristic difference 
between the nPI-method and the standard perturbation theory is that in the lat- 
ter the solutions are written as an expansion in a small parameter, say coupling 
constant, while in the former the equations themselves are expanded. This dif- 
ference is of crucial importance in nonequilibrium quantum field theory; it is the 
very reason for the problems of e.g. secularity with the standard perturbation 
method. The truncation procedure using the nPI-method is practically feasible 
as it turns out that the higher order nPI effective actions become redundant, 
once the order of expansion is fixed. That is, for example in a loop expansioro 
at m-loop order all nPI effective actions with n > m are equivalent. In addition 
to these reductions there may be further simplifications depending on special 
conditions, such as vanishing of the average field [78]. 



3.2.1 From a generating functional to 2PI effective action 

In this work we will concentrate on the two-particle irreducible (2PI) effective 
action [83,88,90-94], which will lead to a self-consistent dynamics for the 2-point 
correlation function Q as well as the 1-point function (0). We show, following 
ref. [9], how the 2PI effective action and the corresponding equations of motion 
are derived for a real scalar field. For fermions we will only give the appropriate 
results. We start by defining 2PI generating functional on the closed time path: 

Z[J,K] = = j Drp[r{0,x)] exp {t{S[r] + Ja^ + ^i^afe^ 0") } , 

(3.10) 

where we use notation with a branch doublet 0" = (0"'", ), so that e.g. Dcf)"" = 
D(j)~^D(j)~, and define a "metric" Cab = cliag(l,— 1), so that Ji(x) = J^(x) 
and J2{x) = — J^(x) and repeated indices are summed over. The CTP ac- 
tion is defined as 5'[0'^] = S[r] — S[(f)~]*, and Ja{x) and Kab{x,x') are lo- 
cal and nonlocal Gaussian sources, respectively. We also use de Witt sum- 
mation convention to leave out integrals in the notation for the source terms: 
Kab^r = J d'^x d'^x' Kab{x,x')(j)°'{x)r{x') and similarly for Ja0". All ra-point 
Green's functions are obtained from Z[J,K] through functional differentiation 
with respect to sources J", while W[J, K] generates the connected n-point func- 
tions. Especially, the average is defined as 



■^A loop expansion in the nPI effective action corresponds to an expansion of the equations, 
not to the perturbative expansion of the .solutions. 

^This is usually called mean field, but in this work we have a different notion for mean 
field, explained later in chapter H) 
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If we set Ja = Kah = after tlie variation, then cj)^ = = (0) is the physical 
expectation value without sources. The 2-point functions can be obtained form 
Vr[J, i^] either through a double derivative with respect to J a or through a 
derivative with respect to the nonlocal source Kab- For later purposes we use 
the latter way, and define the full propagators A'** as 

/^^/'^!, = \ [r + A'^''(x, x')] . (3.12) 

dKab{x,x') 2^ 

From the definition we see that A"^ = (Tc [(0" — 0")((/)'^ — 0")]) i.e. it corre- 
sponds to fluctuations with respect to the average field, and thus it actually 
reduces to Eq. (13.41) only for vanishing average field (f) (and vanishing sources). 
To proceed, we define the 2PI effective action as a double Legendre transforma- 
tion oiW[J,K]: 

r2Pi[0, A] = W[J, K] - J J- - ^Kab [r^' + A-'] , (3.13) 

where it is understood that the sources Ja{x) and Kab are eliminated through the 
relations between them and the correlators 0" and A"^ arising from Eqs. (13.111) - 
(13.121) . These relations are always invertible following from the general properties 
of Legendre transformation. The desired equations of motion for the correlators 
0" and A"^ are now obtained by functional differentiation: 



5A^ ~ ~2^ 



-Ja-Kab4>\ ^ = - Kab, (3.14) 



so that in the case of physical (sourceless) dynamics we get the equations: 
(5r2Pi/50" = and 5r2Pi/5A"'' = 0, which corresponds to finding the extremum 
for the effective action r2Pi[0, A]. 

3.2.2 Formula for the 2PI effective action in terms of the 
fluctuation field 

In order to use the equations of motion (13.141) . we want to find a practical method 
to compute the 2PI effective action r2Pi. To implement the so called background 
field method, we write the effective action in the form: 

e*r-i = J Dr exp[t[s[r]+Jair-r)+lKab{r<P'-r^'-A''')]} , (3.15) 

which follows directly from Eqs. (I3.10p and (I3.13p . Here we have left out the 
initial density matrix contribution p[0"(O,af)]. The justification for this can 
be seen in two ways. First, if the initial state is Gaussian, the initial density 
matrix can be written as p[0"(O,x)] = exp(— |i?afc0"0^), and the new initial- 
time Kernel Rab can be absorbed in the source Kab- But this inclusion seems 
to ruin the desired condition that the physical evolution is given by vanishing 
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sources Ja and Kab- However, since Rab vanishes for all but initial time, we 
see that it affects only the initial conditions for the 1- and 2-point functions 
0" and A"^, and hence can be neglected in the dynamical equations if we do 
adjust these initial conditions correctly [78]. The other possibility when the 
neglection of p [0*^(0, x)] is justified is to consider such initial conditions, where 
the initial state in the distant past ~oo is in the in vacuum. This condition 
is implemented by just shifting the mass rn? to — ie in the first first branch 
and to + ie in the second branch^ i.e. "tilting" the time path in the complex 
plane in the same way as in the standard vacuum quantum field theory [9]. For 
an initial state that is neither of these cases the omitting of p[</)"(0,a;)] is not 
strictly justified and the following developments in this section provide only an 
approximation. If one wishes to consider those non-Gaussian initial states more 
accurately one needs to use higher nPI effective actions. 

To come back to equation (13.151) . we see that using Eqs. (I3.14p and the 
symmetry of the source Kab the exponent becomes 

sm - ^('^'^ - r) - ^ [ir - - 0') - . (3.16) 

Next we shift the integration variable in Eq. (I3.15P by the average fields: 0" = 
0a _|_ expand the classical action in powers of the new fluctuation field 

+ V"] = S[r] + S^aV"" + ^^,a,fe</^V' + ^2 , (3.17) 

where S^a nieans functional derivatives of S with respect to 0" evaluated at 
0" = 0" and similarly for the second derivative 5a,b; and 5*2 denotes the collec- 
tion of the higher order terms (cubic and so forth) in the fluctuation field ^p"". 
Furthermore, we separate the trivial (lowest orders) and nontrivial parts of r2Pi 
by writing it in the form 

r2Pi[0, A] = S\r] + ^2A-^,(0)A"'' - ^zTr[ln A] + r2[0. A] + const , (3.18) 

where the infinite constant (often discarded since it does not affect the equations 
of motion) is — | / d'^x5{Qi) and we denote iAQ^^(0) = Plugging these 

expressions in Eq. fl3.15p we find that the nontrivial part r2 can be expressed 
as: 

= [det A] -1/2 y j^^a 

X exp I - iA,-,VV + ^[S2m - JaV" - Kab{^'^^' - A'^^)] } , 

(3.19) 

^This is why the complex conjugate is exphcitly written in the second branch action 
even though the classical action is always real. 

^The motivation for this notation is that A(^^j(0) becomes the inverse free propagator in 
the limit of vanishing (ff . 
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where 



J a — n^,a. 



b,c 



V 1 n ' 



K, 



ab 



ST2 



(3.20) 



We see that despite the i^abA"''-term r2 has the form of a generating functional 
for a new theory with classical action |A~^^ + 5*2 [(y?"] and sources Ja and 
Kab- Next we show that these sources with the additional term Kab^"'^ will just 
fix the 1- and 2-point functions of this new theory. To show that, we start by 
considering the matrix 




SKab 



(3.21) 



which we know is invertible, because of the general invertibility of Legendre 
transformations. In terms of the derivatives of r2 this becomes (also subtracting 
a singular matrix which does not affect invertibility) 



5^r2 
2 



■ - 2 

i!r2 



<5r2 



(3.22) 



On the other hand, by taking variational derivatives of Eq. (13.191) we get the 
equations 



2 2 



(5A"''(5A<='* 



(^0 







So, since the coefficient matrix is invertible, we conclude that 



(3.23) 



(3.24) 



This result has tremendous implications: As the sources Ja and Kab just kill 
the 1-point function and fix the 2-point function to A'^'^, it follows that we can 
neglect these sources in the practical calculations, if we include only the vacuum 
contributioiill to the effective action r2 using A^*^ as the full propagator of this (p- 
theory. Furthermore, because the /it// propagator is fixed to A'^'^, it follows that in 
the diagrammatic calculations we need to consider only two-particle irreducible 
(2PI) graphs i.e. graphs that do not become disconnected while cutting two 
internal lines (see Fig. 13. 2p . hence the name for the 2PI action. So, we conclude 
that the 2PI effective action for the original theory is given, besides the terms 
explicit in equation (13.181) (classical and one- loop contribution), by the sum of 
all 2PI vacuum graphs in a theory with action |A~^^(y9"y9'' + 5*2 For example, 
for the real scalar field with quartic interaction: 



S[(P''] = S[^+] - S[4>-] = / d^x 



(3.25) 



''in non- vacuum graphs the external legs are connected to the average field of 93-theory that 
vanishes 
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Figure 3.2: Three-loop vacuum graphs in a theory with quartic interaction. Only 
the left one is two-particle irreducible. 



where him = —^2222 = 1 and the other components are zero, we find by per- 
forming the shift 0" = 0" + if"- that the interaction part for the 93-theory is given 
by 

^2^1 = jd'x[- ^habcdr¥>''^V - ^hatcd^V¥>V] ■ (3-26) 

Note that a cubic interaction with an effective vertex depending on the average 
field 0" is generated. 

3.2.3 Schwinger-Dyson equations for the propagators 

Using the expression (13.181) for r2Pi in the equations of motion (13.141) we get 
the following equation for the propagator A'^* in the (physical) case of vanishing 
sources: 

- \'^^>.y) - k'^-Ji^.y) + = 0. (3.27) 



5A'^^(x,?/) 2 "^'""^ 2 '''^ (5A«^(x,i/) 

where the second term follows from variation: 5Tr[ln A]/5A"^ = A^^^. By defin- 
ing the self energy: 

we see that this equation (I3.27P is of the form of famous Schwinger-Dyson equa- 
tion for the full propagator: 

Kb^^^y) = ^oU^^y) + ^ab{x,y) , (3.29) 

which (upon inverting) is presented graphically in Fig. 13.31 To get this equation 
in a form that is feasible for practical calculations we multiply it from the right 
by A"'' to obtain 

j d^zA^^^^ix, z)A''\z, y) = 5ab5\x - y) + J d'^zUadx, z)A''\z, y) . (3.30) 
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Figure 3.3: A generic form of a Schwinger-Dyson equation for the full 2-point 
function A'*''. 

This is an integro-differential equation for the full propagator A"*, because the 
inverse free propagator Aq ^ on the LHS contains explicit spacetime derivatives 
acting on A"^, while on the other hand the full propagator appears inside the in- 
tegral on the RHS, and also the self energy 11 is (typically) a nontrivial functional 
of A"''. One should stress that this Schwinger-Dyson equation is formally exact 
in the case of Gaussian (or distant past vacuum) initial conditions. However, 
the actual computation of the self energy via the 2PI action r2 is a nontrivial 
task and in most cases of interest necessarily involves some truncation, such as 
the loop expansion or a large N expansion for (9(A^)-invariant theories. 

For fermions the construction of the 2PI effective action is very similar to 
scalar fields. However, due to the fact that fermionic fields are Grassmann 
numbers, there are some differences in the final formula, which for vanishing 
fermionic average field is given by [78]: 

T2Pi[G] = -iG^^^.G'"' + zTr[lnG'] + T^IG] + const , (3.31) 

where r2 is again the sum of two-particle irreducible vacuum graphs in the 
theory with propagator G"'' and the vertices are the same as in the original 
action, because of the vanishing average field. By comparing this with the scalar 
expression (13.181) we see a difference of a factor —1/2 in the first two terms. This 
is due to different values of the functional determinant arising from the scalar 
and fermionic path integrals. In the same way as for scalar fields we obtain a 
completely analogous Schwinger-Dyson equation for the propagator G"''^: 

j d^zG^l^{x, z)G'\z, y) = 5ab6\x - v) + j d^'zT.acix, z)G'\z, y) , (3.32) 

where the self energy S is now defined as: 

Note that for a complex scalar field the 2PI effective action would be similar 
to the fermionic expression (I3.3ip . without 1/2 factors, but with different signs 
in the first two terms (see e.g. [37]). Consequently the definition of the self 
energy would be the same as in Eq. (13.331) except for the sign difference. The 
reason why the complex scalar field differs from the real one is basically just the 
doubling of the degrees of freedom for a complex field. Finally, we note that for 
an interacting theory that combines fermionic and scalar fields e.g. by a Yukawa 
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interaction, we would just need to combine the nontrivial r2-parts of the 2PI 
effective actions and include contribution from the Yukawa interaction. This 
would naturally lead to cross-couplings in the fermionic and scalar Schwinger- 
Dyson equations. 

3.3 KadanofF-Baym transport equations 

Let us use again the generic notation Q = {G, A} for fermionic and scalar 
propagators and denote both self energies simply by S, and also adopt the same 
notation: S++ = etc., as for scalar and fermionic propagators Eqs. (13.41) - 
fl3.5p . It follows that the different components of Schwinger-Dyson equations 
(13.301) and (13.321) are consistent provided that the self energies can be divided in 
local (singular) and nonlocal parts 

^^\u, v) = c^'^S^u - t;)Ssg(M) + S"^(m, v) , (3.34) 

where the nonlocal part obeys similar relations as the propagators: 

T,f{u,v) = 9{uo — vo)T,'' {u,v) ± 9{vo — uo)T,'^{u,v) 

T,p{u,v) = 6{vq — uo)T,'' {u,v) ± 9{uo — Vo)T,'^{u,v) (3.35) 

for scalars and fermions, respectively. These relations (I3.34p - (l3.35p should hold 
generally for any reasonable approximation of the self energy [37]. The singular 
term Egg can be absorbed into the inverse propagator Qq^i, on the LHS of the 
Schwinger-Dyson equations either to the mass renormalization or to a classical 
background field (for example for gauge interactions). From now on we assume 
that this absorption is made and denote simply: E"* = S"^. We further define the 
retarted and advanced self energies and their (anti)hermitian parts analogously 
to the propagators in Eqs. (13.61) and (I3.8jfl The antihermitian part is denoted 
as F: 

F = 1 (S'^ - SO = ^ (S> T S<) , (3.36) 

corresponding to the scattering width of the field excitations. For scalar fields 
Eq. (13.361) is conventionally defined as ujT, when F is directly the scattering 
width with the correct dimension. For fermions F is a 4 x 4 (spinor) matrix and 
the physical meaning of various elements is more obscure and will be discussed 
later in section 15.31 in the case of interaction with a thermal background. 

Next we want to write the Schwinger-Dyson equations in a different form to 
make a separation between the dynamical and spectral properties of the system 
more evident. Using the above definitions and the corresponding ones for the 
propagators in Eqs. (13. 6p and (13.80 . and the fact that for vanishing average fields 
the inverse free propagator obeys GqIi, = CabGop, it is a matter of simple algebra 

*The defined self energies clearly follow the hermiticity properties of the propagators, 
Eq. dSUl) and below 
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to show that the Schwinger-Dyson equations fl3.30p and fl3.32p can be written in 
the form: 

{go}-j:H)(»gH + T(^A = 6 (3.37) 

and 

{g-^p -j:j,)(g)g< -j:< (g)g^ = ^ (s> ® 6;< - s< ® g>) , (3.38) 

where, as stated before, we have assumed that the singular self energy Ssg(M) is 
absorbed into gQp, and we use the notation for the convolution integral: 

f®g = j d^zf{u,z)g{z,v). (3.39) 

The equations fl3.37p are called pole equations, while Eq. (13.381) is one of the 
two Kadanoff-Baym (KB) equations [95]. The similar KB-equation for the other 
Wightman function g^ needs not to be considered, since from the definition (13.81) 
it immediately follows that g^ = ±g^ —2iA. In general, the pole equations will 
fix the spectral properties of the theory, while the KB-equations will give the 
dynamical evolution, i.e. the quantum transport effects. Indeed, in the classical 
limit the KB-equations (13.380 for fermions and scalars will reduce to well known 
quantum Boltzmann transport equations for the phase space number densities 
(see e.g. [9,34,36-38]). 



3.3.1 Mixed representation and gradient expansion 

If there is a clear separation between internal (microscopic) and external (macro- 
scopic) scales in the system, it is appropriate to analyze the pole- and KB- 
equations (I3.37I) - (I3.38I) in so called mixed or Wigner representation, where a 
partial Fourier transformation with respect to the internal coordinate r = u — v 
is performed. This transformation leads to a gradient expansion in derivatives 
of the external (average) coordinate x = [u + f)/2, which contains, in general, 
infinitely many terms. However, if the separation of the scales is manifest, this 
expansion can be truncated (or resummed) to a good approximation. 

To begin with, let us define the Wigner transformation for an arbitrary 2- 
point function F{x,y): 

f(...)./.Ve-F(. + ./2..-./2). (3.40) 

where x = {u + v)/2 is the average coordinate, and k is the internal momentum 
conjugate to the relative coordinate r = u — v. Using this definition it is easy to 
transform the equations (13.371) and (13.381) into the mixed representation to get 

g,^'pA-e-''>{^H}{A}-e-''>{T}{gH} = (3.41) 
^o:Ffo-e-*^{S^}{^^^} + e-^^{r}M} = 1 (3.42) 
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and 

Qo,fG^ - e-''>{i:H}{g^} - e-*0{S<}{^^,} = Ceoii , (3.43) 
where the colhsion term in Eq. (13.431) is given by 

Cecil = -te~'^ m{g<} - {S<}{^}) , (3.44) 

and the 0-operator is the following generalization of the Poisson brackets: 

0{f}{9} = \ [dj ■ dug - dj ■ d,g] . (3.45) 
The differential operators Q^p are related to the inverse propagators: 

A-^ = e -^d^ +ik-d-rrr'e--^^'^-^^ (3.46) 

Go> = ^ + ^ -m^e-^^^"^-^^ - z7W-^^"-^\ (3.47) 

for bosons/fermions respectively, where d"^ means that the derivative is acting 
on the left to m^, niRj. Note that these operators are not the transformed inverse 
propagators, but follow from the identity: e'^'^ [Qq]?} {F} = G^pF. 

Eqs. fl3.4ip - fl3.43l) are the desired quantum transport equations for the 2- 
point correlation functions Q"^, Qh and A. We are primarily interested in solv- 
ing the Kadanoff-Baym equation fl3.43p for the Wightman function but in 
general, also the pole equations fl3.4ip -f !3^42|) have to be considered because of 
the cross-couplings between the equations. One can see that these equations 
indeed contain the derivatives of the masses m and the self-energies S up to 
infinite order, restraining their use in practical applications unless this gradient 
expansion can be truncated (or resummed) in some reasonable way. 

In the standard approach to quantum kinetic theory the conditions of so called 
quasiparticle or on-shell approximation are assumed, including slowly varying 
[i.e. "nearly" translation invariant) background fields and correlators, as well 
as weak interactions [9]. This approach provide a consistent approximation to 
truncate the gradient expansion in KB-equations fl3.43p to leading order, culmi- 
nating in the derivation of the famous quantum Boltzmann equations. However, 
because of the assumptions made, it follows that the correlators whose dynamics 
we are studying, are "close" to local thermal equilibrium throughout the evolu- 
tion. Especially the information on quantum coherence will be irrevocably lost. 
In the next chapter we start to build an extended approximation scheme that 
incorporates the good features of the standard kinetic approach with easy-to- 
use Boltzmannian-type transport equations, yet including the effects of nonlocal 
quantum coherence. 
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3.4 Physical quantities from 2-point correlation 
functions 



We conclude this chapter by writing down some famihar physical observables, 
including the number currents and the energy momentum tensors for fermionic 
and scalar fields, in terms of the 2-point correlation functions and A^. These 
expressions follow directly from the definitions of the correlators in Eqs. fl3.4ti3.5l) 
written in the mixed representation, so we simply list the results here. The 
expectation value of the fermionic number current is given by 



and for comj>/ea§ scalar bosons we have: 



d^A; 



Ti [-ff'iG<{k,x)] 



(3.48) 



(j^(x)) = -z(9^0V-0ta^0> 



The symmetric (Belinfante) energy momentum tensor (see e.g 
is 



2FzA<(A;,x). (3.49) 
for fermions 



(2 
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i-f^k" + Yk^)iG<{k,x) 



(3.50) 



while the bosonic tensor for real scalar field is given by 



d^k 



k^k^ + \d^,d:-\g^'^{e 



m 



+ i9^)]zA<(fc,x),(3.51) 



where g'^'^ is the standard Minkowskian metric with signature (+, — , — , — ). For 
a complex scalar field we get a result with the last row of Eq. (13.511) multiplied 
by two. These relations demonstrate the importance of the 2-point Wightman 
functions in nonequilibrium quantum field theory. Later, in section 14.61 we will 
use these results to express the observables in terms of the on-shell distribution 
functions. 



No conserved Noether current can be defined for a real scalar field 
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Chapter 4 



Quasiparticle picture including 
nonlocal quantum coherence 

4.1 Extended quasiparticle approximation 

In the literature (see e.g. [97]) a quasiparticle approximation (QPA) usually refers 
to a set of approximations leading to a spectral phase space structure for the 2- 
point correlation functions composed of sharp singular shells with definite 

energy momentum dispersion relations, such as the standard free particle mass 
shell. An alternative, stronger definition for QPA [9] requires that the func- 
tional forms of the free theory propagators are preserved in (QPA) interacting 
theory, with only the mass parameters replaced by effective masses. For both 
of these definitions the necessary conditions for the quasiparticle approximation 
to be justified include weak interactions and slowly varying background fields. 
Moreover, the standard treatment of QPA relies on the assumption that system 
be close to thermal equilibrium and consequently nearly translation invariant, 
the famous example being the derivation of the quantum Boltzmann transport 
equation from KB-equations [9]. 

Our extended quasiparticle approximation (eQPA) scheme relinquishes the 
assumption that the system needs to be close to thermal equilibrium. The key 
observation in our approach is that under the (otherwise) same conditions of 
QPA with weak interactions and slowly varying background fields, the phase 
space of the 2-point correlators ^^'^ contains novel and completely different 
singular shell solutions, in addition to the standard (quasi)particle mass-shell 
solutions. These new fco,z = 0-shell solutions are unavoidably absent if we de- 
mand that the system is near thermal equilibrium, hence their lacking in the 
standard quasiparticle treatments. We will examine these new solutions in detail 
in section 1131 where we interpret them as describing the nonlocal quantum coher- 
ence between "opposite" (quasi)particle excitations. After the complete spectral 
structure of the correlators is discovered, we feed it as an ansatz to the dynami- 
cal equations to find out the equations of motion for the corresponding on-shell 
distribution functions /, including the novel coherence shells. In this way we 
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get an extension to the quantum Boltzmann transport equation to include the 
effects of nonlocal quantum coherence. 

First, we examine the necessary conditions for the extended quasiparticle 
approximation to be justified. These conditions include weak interactions, slowly 
varying background fields, and existence of particular spacetime symmetries. If 
some of these conditions are not fulfilled, the spectral approximation for the 
phase space breaks down. At the end of this section we outline the procedure 
of using the dynamical equations to find the desired equations of motion for 
the on-shell functions. The actual derivation of the appropriate equations for 
fermions and scalar bosons is presented in section 14. 4[ 



4.1.1 Weak interactions 

The limit of weak interactions in the context of quasiparticle approximation 
means F — in Eqs. fl3.4ip - fl3.43p . where F is the interaction width defined in 
Eq fl3.36l) . This limit is taken strictly whenever the phase space properties of the 
correlators are studied. However, when studying the dynamical properties, one 
has to include F to at least leading order to get any thermalization effects. This 
is precisely the way how the Boltzmann transport equation is obtained in the 
classical limit. 

Neglecting the terms proportional to F in Eqs. fl3.4ip - fl3.43p . including the 
collision term Ccou (in general S^''* ~ F), leads to the free field equations except 
for the self energies Hh- It is not completely obvious how these self energies 
should be handled, however. For a controlled expansion in the coupling constant, 
all the self energies need to be treated in an equal footing; if we neglect F ~ g"', 
we have to neglect also Ej^ at the same order. So, it would be justified to retain 
Ejj only if it is of lower order in coupling constant than F. This is indeed the case 
for gauge interactions for example, with J^h ~ and F ~ in the lowest order. 
Another approach is to treat Hh and F completely independently, so that even 
if there is no coupling hierarchy, J^h will be retained in the equations when the 
phase space properties of the correlators are examined. The motivation for this 
non-controlled approximation is that retaining J^h will not ruin the spectral 
structure of the correlators; it merely modifies the dispersion relations of the 
excitations, so that the standard free particles become quasiparticles. 



4.1.2 Slowly varying background, mean field limit 

In general, it is not enough to neglect the terms proportional to F to obtain a 
spectral phase space structure for the 2-point correlators. This can be demon- 
strated for a scalar field with nonvanishing constant dtm? = e > 0, while other 
derivatives of the mass m are vanishing. The spatially homogeneous solution for 
the free field (Sjj is also neglected) KB-equation 03.430 is then: 

zA<(fco,|fc|;t)ocAif L__^iZJ , (4.1) 
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where Ai(a;) is the Airy-function. This is not a singular distribution in momen- 
tum for nonvanishing e, but indeed in the hmit e ^ it reduces to iA'^ oc 
S{k'^ — m?{t)). This example illustrates that in order to obtain spectral phase 
space structure one needs to consider slowly varying background fields. More- 
over, to actually get singular shell solutions all derivatives of the background 
field have to be neglected in the equations fl3.4ip - fl3.43p . This approximation 
of including only the zeroth order gradients of the background is called mean 
field (or adiabatic) limit. The resulting KB-equations for the study of the phase 
space properties of the correlators in combined weakly interacting and mean 
field limits are then 

{l/^-'-^^-mn-imn'' -^h)G< = (4.2) 

{k'^ -^d'^ +ik- d -m^ -YiH)A< = (4.3) 

for fermions and scalars, respectively. The corresponding pole equations are 
completely the same with replaced by 1 on the RHS of the equations for Gh- 

4.1.3 Special spacetime symmetries 

It turns out that even the mean field equations f l4.2l) -f lT3l) do not in general 
have spectral solutions for the correlators. We see this by considering a 1 -|- 1 
dimensional free scalar field, for which the solution for Eq. fl4.3p reads (for k^ ^ 0) 
[98]: 

iA< {k, x) = [e{-k'^) + e{P - m^)] [A{k) cos(g ■ x) + B{k) sin(g • x)] , (4.4) 

where A{k) and B{k) are real functions of k determined by the initial conditions, 
and q is defined as: 

q ^ 2\koW{k^-m^)/k^ (^^, 1^ . (4.5) 

This solution is not restricted to spectral form with support only on singu- 
lar shells; instead it potentially has support everywhere inside the mass shell 
k"^ > w? , or outside the light cone k"^ < 0, depending on the unspecified func- 
tions A{k) and B{k). For an arbitrary x-dependence this is in conflict with the 
quasiparticle approximation, which requires that the phase space structure is 
spectral. However, if we demand for example the complete translational invari- 
ance: dxo^xiA^{k,x) = 0, then q must be zero, implying that k"^ — m? = 0. We 
conclude that there can be support only on the mass shell, and the solution 
is indeed spectral. In the same way, if we demand only time translational in- 
variance, dxoA'^{k,x) = 0, we find two different spectral solutions: the former 
mass-shell solution, but also a nonconstant solution in xi with ki = 0. Later, in 
section 14. 3[ we examine the spectral phase space solutions in more detail with 
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a different approach by subjecting the solutions of Eqs. fl4.2p - fl4.3p on certain 
spacetime symmetries in the first place. At least for fermions, this seems to be 
the only reasonable method, because of the complex spinor structure. In section 
I4.3l we will find that the phase space structure of fermionic and scalar correlators 
is indeed spectral for two particular spacetime symmetries of interest: spatial 
homogeneity and static planar symmetry. 

4.1 A Spectral ansatz for dynamical equations 

The next step in our approximation scheme is to insert the spectral solutions 
as an ansatz in the full KB-equations fl3.43p . Now we are interested in the 
dynamics of the spectral solutions, so we will consider only those of the resulting 
component equations that include direct spacetime derivatives. We will resort 
to some approximations also in this step. As stated before, we are assuming 
the limit of weak interactions. However, now we do not want to neglect the 
terms proportional to T completely, since that would lead to collisionless plasma 
dynamics. Instead we will include only the leading order terms in F. To see 
which terms are actually leading order can be somewhat difficult in practice. It 
has been shown that for a scalar field close to thermal equilibrium, the term 
proportional to Qh on the LHS of Eq. fl3.43p is of higher order in F than the 
dominant contribution from the collision term [37], so in this limit it is justified 
to neglect that term. For more general situations, however, it is not evident that 
neglecting this term would be justified by any simple arguments. For fermions 
this hierarchy has not been shown even for systems close to thermal equilibrium, 
as far as we know. Nevertheless, neglecting the term oc Qh may anyway be 
a good first approximation; for example the well-known Boltzmann transport 
equation is derived in this limit. Further investigations on the role of this term 
(in dynamical equations) are definitely needed to make any conclusion on its 
importance. 

The role of the term proportional to Ejy on the LHS of Eq. (13.430 poses 
another question. Clearly, if the mean field part has been included in the study 
of spectral properties in equations (14. 2p - (14. 30 . the same term should be included 
in dynamical equations as well. The gradient correction^ to this term are yet 
another question, and in the sense of a controlled expansion in coupling con- 
stant, they should also be included. However, the role of these higher gradient 
contributions would probably be the same as the mean field contribution to 
Ejj i.e. to affect the spectral properties by "modifying" parameters like masses 
and momenta in the equations. By these arguments, neglecting higher gradient 
terms of T^h when the dynamics of the spectral solutions are studied, seems to 
be quite well justified. 

On the other hand, it is of course possible to include all of these interaction 

^The higher gradient corrections in SCJ-tcrms do not necessarily correspond to gradients 
of the background field, and thus are not on the same footing as the overall mean field ap- 
proximation 
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dependent terms in Eq. fl3.43p . Based on the above discussion, it is appropriate 
to do this by including everything else than the mean field contribution of Ej^ 
formally into the collision term. The KB-equations for fermions and scalars then 
become 

(^-^^,.-m^e-t^^"^-^^-z7W-^^-^^-SH)G< = C?,, . (4.6) 
[k^^-^d^ + tk-d-m'e-'^^^-^^ -YiH)A< = C;tu, (4.7) 
where the "extended" collision terms are defined as 

C-toXX = e-'0{{J:H-^r}{G<} + {J:<}{GH + ^A})-EHG<. (4.8) 
Cin = e-^Oi{UH-^T}{A<} + {U<}{AH + ^A})-nHA<. (4.9) 

The higher order gradients in the collision term are a delicate issue in our ap- 
proximation scheme. Usually, when slowly varying backgrounds are studied and 
the solutions are close to thermal equilibrium, it is justified to neglect consis- 
tently all higher (than zeroth) order derivatives in collision term, since those will 
necessarily correspond to higher order gradients in the background field. How- 
ever, in our scheme the coherence shell solutions are rapidly oscillating even in 
constant backgrounds, and consequently the higher derivatives in the collision 
terms do not necessarily correspond to higher order gradients in the background 
field and thus it is not a j)non justified to neglect them. In practical calculations 
this gradient expansion has to be truncated, however, unless the different order 
gradients of the self energy terms can be resummed in some useful way. Later 
on in chapter El we show that this resummation is indeed possible for a scalar 
field interacting with a thermal background. 

To summarize, our approximation scheme works as follows: First, we find 
out the spectral properties of the 2-point functions by using the weakly interact- 
ing mean field equations (14. 2p - (14. 31) . Then, we substitute the obtained spectral 
solutions as an ansatz in the full interacting equations fl4.6l) - fl4.7p (with the ex- 
tended collision terms Ccou or the standard ones of Eq. fl3.44p ) to find out the 
equations of motion for the on-shell functions in the presence of collisions. The 
crucial difference compared with the standard treatment based on the quasipar- 
ticle approximation is that we do not assume that the system is close to thermal 
equilibrium at any moment, yet we are using the spectral solutions for the cor- 
relators. This apparent paradox will be settled in section 14. 3[ where we find the 
novel coherence shell solutions that have completely different properties than the 
standard (quasi)particle mass-shell solutions. In what follows, we will neglect 
the Tjh and IIj^ terms on the LHS of Eqs. fl4.6l) - fl4.7p . since we are primarily 
interested in the general structure of the phase space and the dynamics of the 
on-shell functions, and not so much in the specific modifications of dispersion 
relations caused by interactions. Thus, while studying the phase space proper- 
ties, the corresponding equations are effectively reduced to the noninteracting 
mean field limit. 
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4.2 Reduction of the spin structure in fermionic 
equations 

Before we enter the study of the phase space shell structure, let us first simphfy 
the fermionic KB-equation fl4.6l) further in the cases of two spacetime symmetries 
of interest: the spatial homogeneity and the static planar symmetry. To begin 
with, we write the equation (14 .6^ for the hermitian Wightman function, defined 

as 

G<{k,x) = iG<ik,x)-f\ (4.10) 
By multiplying both sides of Eq. (14. 6 p by 7° we get 

(fco +'-dt-a-{k- - 7°mo - iiWm,)G< = 7°^d7° , (4.11) 

where we use the notation: 

mo,5{x) = mn,i{x)e-^^^-^^ , (4.12) 

and we have dropped the tilde in the collision term to denote either the extended 
collision term in Eq. (14.91) or the standard one in Eq. (13.441) . 

4.2.1 Spatially homogeneous case, helicity diagonal 
correlator 

In a spatially homogeneous case the spatial gradients of and m in Eq. (14.111) 
vanish, and consequently the helicity operator h = k ■ S = k ■ 7^77^, where 
k = k/\k\, commutes with the differential operator on the LHS of Eq. (14. lip . 
This implies that different helicity projections 

G<^,{k,x) = PhG^{k,x)Ph', (4.13) 

where Ph denotes the helicity projector: 

P^ = ^(l + hh), PnPh' = 5hh', h = ±l, (4.14) 

do not mix in a noninteracting theorjQ, so that helicity is a good quantum 
number i.e. a conserved quantity. It follows that the helicity off-diagonals couple 
to the dynamics of the diagonal part only through the collision term. In this 
work we will not consider these cross- couplings, but use the helicity diagonal 
part of the correlator: G"^ = '^^=±1 ^hh^ as an ansatz for the interacting theory. 
In the Weyl basis, where the gamma-matrices are given by the following direct 

^By noninteracting theory we mean here that the self-energies S"'' and consequently the 
collision term are vanishing. The system still interacts with the varying classical background 
(giving rise to varying mass). 
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product expressions (both and a* are the usual Pauh matrices referring here 
to chiral / spin d.o.f. respectively): 

70 = pi(g)l, a = -p^®a and 7^ = -p^ ® 1 , (4.15) 
the helicity diagonal correlator can be written as: 

G<^^g<(^^]^{l + hk-a), (4.16) 

where are (unknown) hermitian 2x2 matrices (for h = ±1) in chiral indices. 
By taking a helicity diagonal projection of Eq. (14. lip we get then an equation 
for g<: 

[ko + ^dt + h\k\p^ - mop' + msp^) (7< = , (4.17) 

where the 2x2 matrix Ch is the chiral part of the helicity diagonal projection 
of the collision term: 

Ph (7°^d7°) Ph = CH0 ^(1 + hk-a). (4.18) 

Given that g^^ is a hermitian matrix, it is useful to decompose the equation 
(14.171) into hermitian (H) and antihermitian (AH) parts: 

(H): 2k,g< = Hg< + g<H^ + C^ (4.19) 

(AH): tdtg;: = Hgt-g<H^ + C^, (4.20) 

where 

H = -h\^p'^ + m^p^ - rh^p^ (4.21) 

can be interpreted as a local free field Hamiltonian operator, and = Ch±Cl are 
the hermitian and antihermitian parts of Ch- We see that in the noninteracting 
mean field limit with Ch = and mo,5 = fn^j only the AH-equation contains 
an explicit time derivative of g"^, while the H-equation is a purely algebraic 
matrix equation. For this reason the AH-equation is called a "kinetic equation" 
describing the dynamical evolution of the helicity diagonal Wightman function in 
a varying background. The H-equation on the other hand is called a "constraint 
equation" , which will constrain the solutions of the kinetic equation in the 4- 
dimensional phase space and thus is the one to be used in determining the phase 
space properties of the (eQPA) interacting theory. 

4.2.2 Static planar symmetric case, spin- 2; diagonal 
correlator 

In a static planar symmetric case all but one spatial derivative (chosen to be 
z here) of and m vanish. Then, apart from the c5|| • fc||-term, the differen- 
tial operator on the LHS of Eq. (14. lip commutes with the spin in ^-direction: 
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— ^o^3^5_ pg^j^ proceed analogously to the previous section |4. 2. 1^ by first 
performing a Lorentz boost to a frame where this term vanishes [34,36,37]. The 
Dirac representation of the desired boost is found to be 

Sn = sgn(ko)^= . (4.22) 

^j2ko{ko + ko) 

The boosted correlatoiEl 

G<Cko, k,; z) = S\^G<{k- z)S^\ (4.23) 

then obeys an equation 

(fco-«' (A;,-l9,J-7°mo-27V^5 ) G<(fco, k,- z) = S-'^'tCi,,^'S, , (4.24) 

where fco = sgn(A;o)(^o ~ ^f\)^^'^- After this boost the differential operator on the 
LHS of Eq. (14.241) indeed commutes with the spin 5*^, as expected. Analogously 
to the spatially homogeneous case this implies that different spin projections 

G< = PsG<Ps' , (4.25) 

where Pg denotes the spin-z projector: 

Ps = ^il + sS'), P,Ps> = Sss', s = ±l, (4.26) 

do not mix in a noninteracting theory, so that spin- 2; is a good quantum number. 
We again neglect the effects of spin off-diagonals and consider only the spin 
diagonal correlator, which is written in the Weyl basis as: 

G<,^g<^^il + sa'), (4.27) 

where gf are hermitian 2x2 matrices (for s = ±1) in chiral indices. By taking 
the spin- 2; diagonal projection of Eq. fl4.24p we get now the following equation 
for gf: 

(ko + s{k, - ^9,)p3 - mop' + rh5p^)gf = , (4.28) 

where the 2x2 matrix Cg is the chiral part of the spin- 2; diagonal projection of 
the boosted collision term: 

Ps (5|7S°^C,tu7°^ii) Ps = Cs® ^(1 + sa^) . (4.29) 

Now, because the matrix is multiplying the derivative dz in Eq. fl4.28p . the 
straightforward division into hermitian and antihermitian parts does not lead 

■^Whereas transforms conventfonally: S C^S"^, the hermitian correlator = 

iG^^^ obeys a pecuhar transformation law: — > SG'^S. 
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to a convenient separation of the derivatives. However, by first multiplying 
Eq. (14.281) from the left by and only then taking the hermitian and antiher- 
mitian parts gives the desired division: 

(H): -2sKgf = Pgf + gfP^ - {p'C,)+ (4.30) 

(AH): isd^gf = Pg<-g<P^-{p^C,)-, (4.31) 

where 

P = kop' + t(mop^ + m,p^) ^ (4.32) 

and (p^Cs)^ = p^Cs ± Clp^ are the hermitian and antihermitian parts of p^Cs- 
These equations are analogous to Eqs. (l4.19l) -f Hl20l) of the spatially homoge- 
neous case. Again, the AH-equation is called a "kinetic equation" describing 
the dynamical ^-evolution of the spin diagonal Wightman function, while the 
H-equation is called a "constraint equation" determining the phase space prop- 
erties of the correlator. 

4.3 Phase space shell structure 

We now begin to examine the phase space structure of the fermionic and scalar 
Wightman functions iG^ and iA^ in the (extended) quasiparticle limit discussed 
in section 14.11 for the case of S// = Uh = 0. We find out that in addition to 
the standard mass-shell excitations, with the dispersion relation k"^ — |mp, the 
phase space consists of novel singular shell solutions that are located at fco = 
for a spatially homogeneous case and at /c^ = for a static planar symmetric 
case. 

4.3.1 Fermions 

As discussed above, the relevant equations that describe the phase space prop- 
erties of the fermionic Wightman function are the constraint (H) equations 
Eqs. ( 14.191) and fl4.30p . that in the noninteracting mean field limit reduce to 

2^0^?;: = {H,g<}, H = -h\k\p^ + mRp^ -mip^ (4.33) 

-2sk,g< = Pgf+gfP\ P= kop' + i{mRp^ + mip') (4.34) 

for the spatially homogeneous and the static planar symmetric cases, respec- 
tively. 

Spatially homogeneous case 

To further analyze the constraint equation (14.331) for the spatially homogeneous 
case, it is convenient to introduce the so called Bloch-representation for the 
chiral matrix g'^: 

9^ = I {9o + g'p') , (4.35) 
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where are the (chiral) Pauh matrices and g^{k^t) are real functions, because 
of the hermiticity of g^. It is easy to see that in the Bloch-representation 
the constraint equation fl4.33p decomposes into a simple homogeneous matrix 
equation 

Bfg^p = , (4.36) 
where the 4x4 coefficient matrix is (index ordering is here defined as a = 



0,3,1,21 



/ ko h\k\ —rriR mi \ 



B 



(4.37) 



h\k\ fco 
—rriR fco 
\ mi A;o / 

A homogeneous matrix equation, such as Eq. (14.361) . may have nonzero solutions 
only when the determinant of the matrix vanishes. Here the determinant is 
simply: 

det(5) = {e - |m|2) kl , (4.38) 
which implies that the nonzero solutions are possible only when 

A;2-|m|^ = or ko = 0. (4.39) 

These constraints give rise to a singular shell structure for the solutions, since 
they need to be proportional to 6{kQ — A;^ — |mp) or 6{kQ). The former class is 
identified as the standard one particle mass-shell solutions, with the dispersion 
relation ko = ±uJk = ±(fc^ + |mp)^/^, while the latter class of solutions with 
ko = are completely novel in the context of quantum field theory (see Fig. l4.ip . 
Based on the observation that the quantum interference of the plane waves 
~ g±«<^t contains a contribution with ko = 0, we make an interpretation that 
these additional (ko = )-shell solutions describe the quantum coherence between 
the particles and antiparticles (positive and negative energy states) with opposite 
momenta and spin. 

The explicit matrix structure of these solutions is found easily by setting 
ko ^ and /cq = in the matrix equation (14.361) for the mass-shell and the 
coherence shell solutions, respectively [I]. The full chiral matrix corresponding 
to the mass-shell solution is given by: 

1^1; ^ 2./i„(i^u)i.„i ( 1 -^^ ^ ;;/*»^^^ ) - . 

(4.40) 

where Skn = sgn(/co) and (|^|, ^) are real functions parametrizing this solution. 
The on-shell distribution functions are called the phase space densities for 
positive and negative energy modes, respectively. Indeed, in the thermal limit 
they are related to the number densities of physical particles, as we will see later. 
For ko = 0-shell solution we get on the other hand 
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Figure 4.1: Shell structure of the correlators and A*^ in the spatially ho- 
mogeneous case. The dark filled (red) blobs show the mass-shell contributions 
for a given \k\, and the light (yellow) blob shows the corresponding coherence 
contribution from the new ko = 0-shell. 



hmfj/\k\ 1 

1 —hrriji/lkl 

—hmi/\k\ —i 

i hmi/\k\ 



6{ko) , (4.41) 



where fi{\k\,t) and /2^(|fc|, t) are new undetermined real functions corresponding 
to the degrees of freedom of this coherence solution. The most general solution 
satisfying the constraint equation fl4.33p (or equivalently the matrix equation 
(14.361) ) for a spatially homogeneous case is the linear combination of Eqs. (14.401) 
and fICTD : 

(4.42) 



9h 



9h,r 



+ 9h,o- 



This general solution contains four independent on-shell distribution functions 
/± 1 2 (fo^ both h = ±1), which is just the number of independent components 
in a hermitian 2x2 matrix, such as the chiral matrix g^. Indeed, in section 1^7^ 
we find that there is a one-to-one mapping between these on-shell functions and 
the components of the fco-integrated chiral matrix (g^)- 



Static planar symmetric case 

The analysis of the static planar symmetric case proceeds in complete analogy. 
By introducing a Bloch representation for gf: 



1 



9f {90 + 9!p') 



(4.43) 
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Figure 4.2: Shell structure of the correlators and in the static planar 
symmetric case. The filled blobs show the contributions from the mass-shell 
and the new kz = 0-shell for a given k^ = uj. Peaked functions on each axis 
represent the accuracy of possible exterior knowledge or "measurement" on fco 
and kz allowing coherent mixing (wide solid peak in /c^-axis) or reducing the 
system to a noncoherent evolution (dashed peak). 



with real g^, the constraint equation 04.341) decomposes again into a homoge- 
neous matrix equation 

Bj^gl = , (4.44) 
where the 4x4 coefficient matrix is now (index ordering is again a = 0, 3, 1, 2): 



B 



(4.45) 



/ ko skz —rriR mj \ 

skz ko 

rriR skz 

y —mi skz j 

The nonzero solutions of this matrix equation are found at the zeros of the 
determinant 

det(5) = A;2(A;2 - Imp), (4.46) 

which are now: 

fc2_|m|2 = or kz = 0. (4.47) 

The former condition gives again the standard one particle mass-shell solutions, 
with kz = ±fcm = =t(^o ~ |?72,p)^/^, while the latter condition kz = gives now 
a different novel class of solutions (see Fig. 14. 2p . By the same argument as in 
the spatially homogeneous case, we interpret that these additional (kz = )-shell 
solutions describe the quantum coherence between the states of same spin and 
energy travelling in opposite z- directions . 

The matrix structure of these solutions is found by setting kz ^ (mass-shell) 
and kz = (coherence shell) in the matrix equation f l4.44p . For the mass-shell 
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solution we get 

(4.48) 

where = sgn^k^) and /^^^ (/cq, 2;) are real on-shell distribution functions 
parametrizing this solution. The coherence shell solution is given by 



where /f (/^o, z) and fHko, z) are new real distribution functions. The most gen- 
eral solution satisfying the constraint equation f l4.34p for a static planar sym- 
metric case is now the linear combination of Eqs. fl4.48p and (14.491) : 

= 9f,m-s + 9s,0-s ■ (4-50) 

In this case also, we later find that the four independent on-shell functions f±i2 
will uniquely correspond to the components of the integrated chiral matrix (gf), 
where the integration is now over instead of /cq- 

We later see that even in the simplest possible example of constant mass 
and no interactions, the new /co,z-shell solutions for the spatially homogeneous 
and static planar symmetric cases are not constant, but oscillate rapidly with 
the frequencies 2ujj^ and 2km, respectively. Thus they break the translational 
invariance of the correlator badly even in this trivial limit. This is the 
very reason why these coherence solutions have not been found (used) in the 
standard treatments of quasiparticle approximation, where it is assumed that 
the correlator is close to thermal equilibrium. 

4.3.2 Scalar bosons 

Next, we discuss the singular shell structure for scalar fields. It appears that 
the method of finding these solutions is not so straightforward, since there is no 
purely algebraic equation in this case. To analyze the phase space properties of 
the scalar Wightman function iA"^, we consider the KB-equation fl3.43p in the 
noninteracting mean field limit: 

(jt^ -^d^ + ik-d-m^yA< = 0. (4.51) 

Like for fermions, we first consider the spatially homogeneous case and then the 
static planar symmetric case. 



rriR/ko 1 
1 uir/ ko 

—mj/kQ —i 
i —rrii/kQ 



Kkz) , (4.49) 
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Spatially homogeneous case 

In the spatially homogeneous case the spatial gradients of the correlator A*- and 
the mass m vanish. Splitting the equation (14.511) into real and imaginary parts 
gives then 

(k^ -Tn" -^d^yA<{k,t) = (4.52) 
kodtiA<{k,t) = 0. (4.53) 

In contrast to the fermionic case we cannot divide these equations into "kinetic" 
and "constraint" equations, since both of them contain time derivatives. How- 
ever, we can now determine the phase space structure indirectly by using both 
of the equations (14.521) and (I4.53P in an appropriate way. 

We begin by setting k^ ^ 0. Then Eq. (14.531) requires that dtiA^ = at all 
times implying that also dfiA^ = 0. Substituting this to Eq. (I4.52p now leads 
to an algebraic equation 

{k'-m')iA<_, = 0, (4.54) 

which has the spectral solution: 

tA<_^{ko, \k\,t) = 27rsgn{ko)fs,^{\k\,t)6{e - m^) , (4.55) 

corresponding to the standard one particle mass-shell solution, with the disper- 
sion relation /cq = = ±(fc^ + m^)^/^. Note that this solution satisfies the 
equation (14.531) only in the mean field limit i.e. neglecting all the time derivatives 
of the mass m. This is expected, since based on the discussion in section 14. 1[ 
the solutions will spread in phase space when the gradients are taken in account 
(see Eq. SI]). 

On the other hand, by setting fco = in the first place, Eq. (I4.53P is identically 
satisfied, and no extra constraint for the time derivatives of iA"^ follows. Then 
Eq. (I4.52P simply becomes 

(d^ + 4^f) A<_, = , (4.56) 
which has the (mean field) solution, 

zA<_,(fco, \k\,t) = 2n[A{\k\,t)cos{2u;^t) + B{\k\,t)sm{2u;^t)] 6{ko) (4.57) 

where yl(|A;|,t) and B{\k\,t) are real functions that become constants when the 
mass m (and ujj:) is a constant. The 5(fco)-factor explicitly fixes the restriction 
to the shell ko = 0. From now on, we will call (parametrize) the factor in the 
square brackets in Eq. (I4.57P as /c(| ^1,^)0, so that the solution (14.570 is written 
simply as 

iA<_Ako,\k\,t)=2nU\k\,mko). (4.58) 
"^Note that fc is not dimensionless here, but has the dimension of 1/M 
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The full spectral solution satisfying the (mean field) equations fl4.52p and fl4.53l) 
is then the combination of Eqs. f l4.55p and (14.581) (see Fig. 14.11) : 

zA< = zA<_, + tA<_s . (4.59) 

This complete solution has three independent on-shell distribution functions 
f±^c, which are now one-to-one related to the three lowest /co-moments of iA^, 
as we will see in section 14. 4[ 

Static planar symmetric case 

For the static planar symmetric case with dt^x,yiA'^ = dt^x,y'^ = the real and 
imaginary parts of Eq. (14.511) are 

(^k^ -m^ + ^dfjiA<{k,z) = (4.60) 
k,dJA<{k,z) = 0. (4.61) 

The analysis proceeds now in complete analogy with the spatially homogeneous 
case. For fc^, 7^ we have d^iA^ = d'^iA^ = 0, so that we find the same 
mass-shell solution as before: 

tA<_^{ko, \k^^\,k,,z) = 27rsgn(A;o)/s,,(A:o, |4U)5(P -m^) . (4.62) 

A new solution is found by setting fc^ = in the first place, leading to (the other 
equation is identically satisfied) 

{dl + Akliz)) tA<_, = , (4.63) 

where km = (k^ — kf^ — m?Y^'^. Analogously to the spatially homogeneous case, 
this equation has a (mean field) solution that can be parametrized as 

iA<_^{ko,\hlk., z) = 27r f,{ko,\hl z)5ik,) , (4.64) 

where fc{ko, z) is a real on-shell distribution function corresponding to this 
coherence solution. The most general solution satisfying the (mean field) equa- 
tions (14.601) and f l4.6ip is once again the combination of Eqs. f l4.62p and f l4.64p 
(see Fig. USD: 

zA< = A1i-s + ^Ao<-s • (4.65) 
In this case the three independent on-shell distribution functions f±^c are one- 
to-one related to the three lowest fc^-moments of iA^. 

Like for fermions, we interpret the new fco,2 = 0-shell solutions for scalar fields 
as describing the nonlocal quantum coherence between the states with opposite 
momenta and 2;-momenta in the spatially homogeneous and the static planar 
symmetric cases, respectively. For a complex scalar field the ko = 0-coherence is 
between particles and antiparticles, while for a real scalar field it is between the 
different modes of the same particlqj- The oscillatory behaviour of the coherence 
shell solution with the frequency 2a;^ in the constant mass limit is now directly 
seen from Eq. (14.570 . 

^For a real scalar field particles are their own antiparticles 
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4.4 Equations of motion with collisions 



Having determined the singular phase space shell structure in the (extended) 
quasiparticle limit, we now want to use the corresponding spectral correlators 
as an ansatz for the dynamical equations including the gradients and collisions. 
Our goal is to find a closed set of equations of motion for the on-shell functions / 
(or some generic quantities related to these). We will only consider the spatially 
homogeneous case here, but the static planar symmetric case would be highly 
analogous to thi^. To proceed, the singular shell structure of the correlator 
naturally suggests to integrate the dynamical equations over the momentum. 
This procedure is supported also in the sense that (local) physical quantities, 
like currents, energy density and pressure, are obtained from the correlator Q"^ by 
a full 4-momentum integration (see sect ion [37il) . There is one profound difference 
between the procedures for fermions and scalar bosons, however. For fermions, 
it is enough to consider the bare integral (zeroth moment) of the correlator to 
find out a closed set of equations for the on-shell functions fa- For scalars, on 
the other hand, the different moment integrals couple to each other, and due 
to our decomposition with three on-shell functions /±,c we need to consider the 
three lowest moment integrals to get a closure. 

4.4.1 Integrated matrix equations for fermions 

To derive the equations of motion for the on-shell functions /±,i,2 for fermions 
in the spatially homogeneous case, we start with the kinetic (AH) equation 
Eq. ( l4.2Up for the chiral 2x2 matrix g^: 

idtg< = Hg< - g<H^ + , (4.66) 

where the operator H is given by Eq. fl4.2ip and Cj^ is the antihermitian chiral 
part of the collision term: Ph (^7°^C^ii7°j Ph, defined in Eq. (14.181) . Integration 
over fco gives now the following equation: 

idtpu = [H,ph] + {C^), (4.67) 

where we use notation (...)= J ^(. . .) for the /co-integral, ph = {g^) denotes 
the (zeroth moment) integral of the correlator g^, and the local Hamiltonian 
operator is given by 

H = -h\k\p^ + mnp^ - mip"^ . (4.68) 

The /co-derivatives appearing in Eq. (14.211) have disappeared as total derivatives 
due to the integration. When we substitute the spectral correlator Eq. (14.421) 

^For the free theory the static planar symmetric case wiU be considered in connection with 
the Klein problem in section [5TT1 
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as an ansatz for g^, the (chiral) components of the integrated correlator are 
simply related to the on-shell functions fa, « = ±,1,2: 



Ph,LL — 
Ph,RR = 




m* 1 

PH,RL = - f-) + ^UU ^f') ■ (4.69) 

It is easy to show that these linear relations can be inverted whenever \k\ ^ 
0. Moreover, because of the singular shell structure of the spectral correlator 
Eq. (14.421) . the collision integral (C^) is composed of projections of the self en- 
ergies S^'^ on these singular shells, and is thus depending on the correlator g'^ 
only through the on-shell functions /„. Based on these observations, we con- 
clude that the matrix equation (14.671) can be used to find the desired closed set of 
equations of motion for the on-shell functions /„. We will not give the resulting 
equations explicitly here, since they would be more messy than Eq. (I4.67p . and 
in any case the specifications of the self energies in the collision term are needed 
to get the equations in the final form. Instead, we refer to the matrix equation 
(14.671) together with the relations (14.691) as our master equations to study the 
dynamics of either the on-shell functions fa or the components of the chiral 
matrix p^. When the self energy functional S*^'^ are specified, these equations 
can be used in practical calculations. In chapter O we use them to study co- 
herent particle production in an oscillating background field in the presence of 
decohering collisions. 

4.4.2 Moment equations for scalar bosons 

For scalar bosons the relevant dynamical equation is Eq. (I3.43p . which in the 
spatially homogeneous case reduces to the following coupled equations (real and 
imaginary parts): 

{k'-\d^-m'cos{ldrdtjyA< = C+ (4.70) 
(fco9i + m2sin(ia,X))zA< = C" , (4.71) 

where are the real and imaginary parts of iC^^n, respectively. We now see 
that because of the explicit fco-factors in Eqs. (I4.70p - (l4.7ip the fco-integration of 
these equations will couple different moments 

p4\k\,t) ^ {k^zA<) = 1^ k^tA<{ko,\k\,t) (4.72) 
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in the same equation. That is, by taking the zeroth moment of Eq. (14.701) and 
the zeroth and first moments of Eq. (14.711) we find the following coupled moment 
equations: 



dtpi = {C-) 
dtP2-ldt{m^)po = {koC-). (4.73) 

Moreover, when the spectral correlator Eq. (I4.65P is used as an ansatz for A"^, 
the moments Po,i,2 are related to the on-shell functions f±^c in a simple way: 

Po = 7r—if+-f-) + fc 
Pi = lif+ + f-) 

P2 = (4.74) 

We see that this linear set of equations is invertible whenever \k\ ^ 0, and also 
the collision integrals (C^) and {koC~) are depending on the correlator only 
through the on-shell functions f±^c, in complete analogy with fermions. So, we 
conclude now that the moment equations (14.731) can be used to obtain the desired 
closed set of equations of motion for the on-shell functions fj-^c- Again, we will 
not give these equations here explicitly, but refer to Eqs. (I4.73I) -( 1T741) as our 
master equations for the study of the dynamics of either the on-shell functions 
/_!_ ,, or the three lowest moments po,i,2- In chapter [5] we use these equations to 
study the coherent production of unstable particles in an oscillating background. 



4.5 Spectral function and thermal limit 

Let us next consider the fermionic and scalar spectral functions A in terms of our 
eQPA scheme. In our approach these are needed in the evaluation of the collision 
terms in Eqs. (14.671) and (14. 73 p . We consider here only the spatially homogeneous 
case, since all the problems we are going to study with collisions in chapter [5] 
have this particular symmetry. The relevant equations for the study of the phase 
space properties of spectral functions (in the noninteracting mean field limit) are 
the same as the corresponding ones for the Wightman functions and A*~, 
Eqs. (14. 20 - 04. 3p . Consequently, the spectral functions have spectral solutions 
given by Eqs. (I4.42p and (I4.59p . with yet undefined on-shell functions f±\2 
fermions and f^^ for scalars. In addition, however, the spectral functions must 
obey the sum rules that can be derived from the pole equations (I3.41l) -( l3^^2ll 
using the spectral relation (13. 9p (or alternatively from the canonical equal time 
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(anti) commutation relations of the fields): 

j ^^(A;,x)7° = l (4.75) 

for fermions, and 

-^{k, + -dt)A{k,x) = l (4.76) 

for scalars. It is easy to show that for the fermionic spectral function the con- 
ditions from the sum rule (14.751) completely fix the values of all the on-shell 
functions 

fl-^ = \ and /f/ = (4.77) 

for both helicities, so that the spectral function reduces to the familiar local 
thermal equilibrium form: 

A = TTSgn{ko){p + mR- i'y^mj)6{P - |mp) . (4.78) 

For the scalar fields, on the other hand, the sum rule fl4.75l) does not immediately 
fix the values of the on-shell functions f^^. However, by use of the dynamical 
equations for the moment functions Pn= J ^ A that are identical to fl4.73l) 
with vanishing collision terms, the sum rule gives that P'o = P'2 = ^ and = 1/2 
in order to get a continuous constant mass limit. The connection relations 
identical to Eq. f l4.74p between p;^ and f^^ then give: 

fi = l, // = 0, (4.79) 

so that the spectral function A again reduces to its standard thermal form 

A = 7rsgn(A;o) 5 - m^) . (4.80) 

Note that the contributions from the coherence shell ko = are completely 
absent in the (eQPA) spectral functions (I4.78P and (14.801) . because of the strong 
constraints imposed by the spectral sum rule. 

Let us conclude this section by giving the Wightman functions and in 
the thermal limit. It is easy to show that the full translation invariance in ther- 
mal equilibrium kills the coherence shell contributions completely, and further- 
more implies the Kubo-Martin-Schwinger (KMS) conditions for the Wightman 
functions [99-101]: 

g>^(ko) = eP''G<^{ko) . (4.81) 

The KMS-conditions and the relation = — 2iA as well as the thermal 
spectral functions Eqs. (14.781) and (I4.80p will then fix the mass-shell distribu- 
tion functions to Fermi-Dirac and Bose-Einstein distributions for fermions and 
bosons, respectively: 

/io(i^i) 



^fq(^o) 



e/3fco + 1 



(4.82) 
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leading to the following thermal equilibrium correlators [97]: 

zG< = 27isgn{ko){p + mii-i'y^mj)n^^{ko)5{k'^ -\m\'^) 

2A< = 27rsgn(A;oX(M'^(^'-"^') • (4-83) 

These expressions will be used in chapter E] when we compute the collision terms 
in the case of interaction with a thermal background. 

4.6 Physical quantities in terms of the on-shell 
functions 

Having found the spectral phase space structure of the correlators and A*~, 
we want to write some of the physical observables considered in section 13.41 in 
terms of the on-shell functions /. In addition, we introduce the concept of phase 
space particle number. 

4.6.1 Particle number and fluxes 

In contrast to standard vacuum QFT, the concept of particles is not very well 
defined in nonequilibrium quantum field theory [8,9]. The problem is that dif- 
ferent observers will see different vacuum states, and on the other hand vacuum 
states are not remaining "empty" even in the free field evolution, because of 
the nontrivial background. That said, in the case of clear asymptotic regions 
in the background field, the problem of particle production from the classical 
background is well understood. Most definitions for particle number are rely- 
ing on some sort of diagonalization of the Hamiltonian [8,9, 102, 103], the idea 
being that the energy of the k-mode would entirely consists of a vacuum part 
plus (anti)particle excitations: Sj: oc oj^{n^ + =F 1) for fermions and bosons, 
respectively. This idea is of course supported by thermal equilibrium, where this 
relation is satisfied. 

Here we adopt a somewhat different concept for the particle number. At the 
end of the previous section we saw that in thermal equilibrium the mass-shell 
distribution functions of and A*^ are just the Fermi-Dirac and Bose-Einstein 
distributions, which correspond to particle number densities for positive energies: 
n^'^ = n^q^{ko = Uj:), and to antiparticle number densities for negative energies 

according to the Feynman-Stuckelberg interpretation: n'^'^ = ±1 — n^^(/co = 
—ujj^)- Within our eQPA scheme the Wightman functions G'^ and A^ always 
have a spectral phase space structure including the singular mass-shell, even in 
out-of-equilibrium conditions. This suggests us to generalize these relations and 
use them as definitions for the phase space particle and antiparticle numbers, 
i.e. for a spatially homogeneous system we define: 

r^kkit) = f1{\k\,t) and n^^^{t)^l-fi{\k\,t) 
n^{t) ^ Ui\k\,t) and ^ _l _ /_(|^|, t) (4.84) 
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for fermions and scalars, respectively. The soundness of these definitions is 
consohdated by expressing the number current densities in Eqs. (13. 481) - ( 13^491) in 
terms of the on-shell functions. We get 

J^siU + f-) = j - H - 1) ' (4-85) 

which are just the expected results with the correct vacuum energy contributions. 
Now, using the inverse relations of Eq. (I4.69P we can express the fermionic 
(anti)particle numbers in terms of the Bloch components of the integrated chiral 
matrix ph = |((^o) + (9^) " 

n^uh = ^{'h\k\{g's)+rnn{g',)-mr{g!^)) + ^{g^,) 

= ^[~h\k\{gi;)+mM)-mj{g!^))-^{gl;) + l. (4.86) 

For scalar fields, on the other hand, the (anti)particle numbers can be related 
to the three lowest moments Po,i,2 by the inverse relations of Eq. (14.741) : 

1 

n,^ = —P2 + P1 

= — P2-P1-1. (4.87) 

We see that the particle and antiparticle numbers coincide if {g^) = 1 for 
fermions and pi = —1/2 for scalar^. These values thus correspond to zero 
"chemical potential", for which the number current densities (14.851) have only 
vacuum contributions. 

Let us now compare our phase space particle number to other definitions 
in literature. Apart from some sign conventions our fermionic particle number 
density in Eq. (I4.86P with zero chemical potential reduces to the one obtained in 
ref. [102] by using the Bogoliubov transformation to diagonalize the Hamiltonian. 
Moreover, using the free field equations of motion (14.731) . we see that our scalar 
particle number in Eq. (14.871) agrees with the one in ref. [103], namely 

{^k + = P<^ H-^o + Id^Po) , (4.88) 

if we take the adiabatic limit: dfpo <^ ^|Po- These comparisons further strengthen 
the soundness of our definitions. 

^For a real scalar field particle is its own antiparticle, thus pi — —1/2 is the only sensible 
value 
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For spatially dependent problems one is often more interested in fluxes, 
i.e. the spatial components of the number currents, than densities. Using 
Eqs. ( I3.48p -( l3ll9l) we find that the z-fiuxes in a static planar symmetric case 
are given by: 

(jUz)) = /7||^^g-(M(/+-/-)- (4.89) 

In this case it is natural to interpret or define f± as right/left moving parti- 
cle fluxes per unit volume in (/cq, fc||)-phase space. We use these definitions in 
sections 15.11 and 15.21 for solving quantum reflection problems with our methods. 



4.6.2 Energy density and pressure for spatially 
homogeneous systems 

Let us proceed by considering the energy density and pressure in a spatially 
homogeneous case. These are by definition the 00- and ii-components of the 
energy momentum tensor, respectively. Using Eqs. fl3.50p - fl3.5ip we get for the 
energy density (in terms of n and n instead of f±): 



d^k 1 



-UJ] 



(27r)3 2 

for fermions and scalars, respectively. For the pressure we get instead: 



(4.90) 



{P^t)) = {d%t)) 
{Ps{t)) = {T\t)) 



d^A; 1 



(27r)3 3 [u^. 



k 



— i^kh + ^kh - 1) - ^Rfi + ^if: 



d^k 1 



(27r)3 3 



Pi 
^2 



K^ + % + i)-(3^I-^')/c 



. (4.91) 



The expressions for energy densities in Eq. fl4.90p expected, consisting 

of particle and antiparticle contributions and the vacuum energies, which are 
opposite for fermions and scalars. The expressions for pressure in Eq. f l4.9ip . 
however, involve explicit contributions from the coherence shell as well. Thus, in 
the presence of coherence the "quantum pressure" may be drastically different 
from the standard classical one. However, because of the typical oscillatory 
behavior of the coherence solutions at microscopic time-scales Atosc ~ l/w, this 
deviation would average out in most cases of interest at time-scales much larger 
than Atosr- 
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4.7 On the validity of eQPA approximation 



We conclude this chapter by discussing the main concerns regarding the vahdity 
of our eQPA approximation scheme, the breakdown (spreading) of the singular 
shell picture and the handling of the gradient expansion in the collision term. To 
begin with, we recall that the basic requirements of the quasiparticle approxima- 
tion, weak (quantum) interactions and a slowly varying background field, should 
be fulfilled for the singular shell picture to be a good approximation. However, 
in our eQPA scheme we relax the assumption that the system is close to a local 
thermal equilibrium, and hence nearly translation invariant. This allows us to 
include the solutions describing nonlocal quantum coherence, which are oscilla- 
tory (and thus translationally non-invariant) on quantum scales with frequencies 
of order 2uj^ (or 2km in a static planar symmetric case). To that end, it is worth 
noticing that in this work we have considered the coherence only in the cases 
of two particular spacetime symmetries: (/cq = 0-shell) in a spatially homoge- 
neous case and [k^ = 0-shell) in a static (or stationary) planar symmetric case. 
For other spacetime symmetries, it is expected that similar coherence solutions 
would exist, living in different singular shells in the phase space. One interesting 
symmetry, worth studying in the future, is the rotationally symmetric geometry. 

Of course, the singular shell picture is not exact for any (nontrivial) practical 
application, as the shells will spread out in the phase space when the interaction 
width r and the gradients of the background field are nonvanishing. Thus, an im- 
portant question concerning the eQPA scheme is whether the new /cq.z = 0-shells 
are spreading similarly (in magnitude) to the well-known mass-shell solutions, 
when the strict quasiparticle limit is lifted. To elaborate on this question, we 
review a simple nonrelativistic example from ref. [Ill], where we computed the 
exact free-field Wightman function from the one-particle wave functions in the 
case of an infinite square well potential. In this case, the standing waves of bound 
states are highly coherent superpositions of the opposite travelling plane waves. 
The essential phase space structure is plotted in Fig. 14.31 We clearly see that the 
spreads of the mass-shell contributions (at k^ = ±A;„) and the k^ = 0-shell are 
of the same order. This observation is confirmed by an analytical calculation, 
which shows that the functional forms of individual peaks are actually the same 
in this specific case. That is, the plotted phase space density in the middle of 
the square well is given by 



where the density corresponding to each "shell" is f{q) = sin(gL)/g. Although 
this is just one simple example, it supports the expectation that the different 
phase space shells are spreading in an equal manner also more generally. In- 
deed, the spreads for all shells should be controlled by the same parameters, the 
interaction width F and the gradients of the background. 

Another important issue in our eQPA scheme is the treatment of the gradient 




(4.92) 



± 
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Figure 4.3: The phase space structure of the 2-point Wightman iA^ in the 
case of a 1-dimensional infinite square well potential. Shown is the phase space 
density Eq. f l4.92p as a function of k^/kn at the middle of the square well at 
z = 0, for odd bound states with n = 9 (left) and n = 99 (right). For more 
details, see ref. [III]. 

expansion in the collision term. For the flow term this expansion is trivial as the 
/c-derivatives can be handled easily in the /c-integration (partial integrations are 
straightforward). For the collision term, on the other hand, nontrivial contribu- 
tions arise that are proportional to spacetime derivatives of the correlator d!^Q^. 
In the standard kinetic approach these terms would necessarily be proportional 
to r or the gradients of the background field c}"?7i[^, as the correlator is adia- 
batic. In our eQPA approach, in contrast, this is not the case: because of the 
oscillatory coherence solutions, one typically has (in a spatially homogeneous 
case): dfQ"^ ~ —Auj'iQ^, such that all orders of the expansion involve terms that 
are not controlled by any small parameter. This severe-looking problem can be 
cured, however, by a recursive use of the equations of motion and a resummation 
of these coherence contributions, to get a controlled expansion in powers of F 
and d'^m. In an example of coherent production of decaying scalar particles in 
section 15.41 this expansion is worked out for a thermal interaction to the leading 
order, neglecting terms of order 0{r'^,r dm) and higher. 



^This notation represents the n-th gradient order, where it is assumed that the n-th deriva- 
tive and the first derivative to the power n are of the same gradient order 



50 



Chapter 5 
Applications 



5.1 Klein problem 



As our first application for tlie eQPA sclieme introduced in cliapterlHwe consider 
tlie famous Klein problem i.e. quantum reflection from a step potential. The 
problem with noninteracting quantum fields can be solved conveniently using 
the Dirac and Klein-Gordon equations for fermions and scalars, respectively. 
However, here we want to illustrate the necessity to include the coherence shell 
solutions when the problem is solved with the methods of quantum transport 
theory. Moreover, when the interacting fields are considered, the wave equation 
approach is not possible anymore, and one needs to use some more advanced 
methods, like our approximation scheme, to solve the problem. 

The setup for the problem is illustrated in figure 15.11 We consider a step 
potential 



such that the mass-shell momenta are k = k]^ = {kl — |mp)^/^ and q = k^ = 
((/cq ~ V)'^ — |mp)^/^ in regions I and II, respectively. The momentum k is 
always real, while q can be either real or imaginary, depending on the values 
of energy ko and potential V. The latter case describes the quantum tunneling 
inside the potential barrier. Without collisions energy is a conserved quantity, 
and we fix now: kQ = lv > 0. In our approach the particle fluxes are just the 
mass-shell distribution functions f±. The asymptotic boundary conditions for 
this setup are the following: We normalize the incoming flux from the right to 
unity /i = 1, and we have no incoming flux form the left = 0. Also, since 
there is no incoming flux from the left, we set the coherence asymptotically to 
zero in region II. In addition, in the case of imaginary mass-shell momentum q 
in region II, we cannot have any asymptotic mass-shell solution there, since no 
propagating wave would penetrate a potential barrier to infinite distance. The 
reflected and transmitted fluxes we want to compute are and respectively. 
We will first consider fermions and then scalar bosons. For both, we analyze the 
cases with real and imaginary q separately. 




(5.1) 
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Figure 5.1: Reflection from a step potential Vi^z) = V6{—z) with the momenta 
of incoming (red) and outgoing (green) particles described by arrows. 

5.1.1 Fermions 

For a varying potential instead of the mass, the fermionic free field KB-equation 
becomes in the mixed representation 

( ^ +1 k - ^^^^i - loV{x)e-^^^-^'')G<{k, x) = . (5.2) 

In this case the potential is planar symmetric, and the reduction of the spin 
structure is completely analogous to section 14.2.21 leading to the coUisionless 
equations 

(H): -2sKgf = Pygf + gfPl (5.3) 
(AH): isd^gf = Pyg<-g<Pl (5.4) 

with 

py = {ko-V)p^ + i{mnP^ + mip^), (5.5) 

where we denote V{z) = V{z)e^^^^'"' . The spectral phase space structure of the 
correlator gf that is obtained by analyzing the H-equation (15.31) in the mean field 
limit, is now identical to the case of a varying mass given by Eq. f l4.50p . except 
for the replacement ko —>■ ko — V. By integrating the kinetic (AH) equation (15. 4p 
over kz we get 

isd^ps = Pvps - PsPy > (5-6) 

where we denote Ps = J ^9f- The fc^-derivatives in V have disappeared due 
to the integration reducing V and Py to their corresponding mean field forms 
V and Py. When the spectral correlator is used as an ansatz for gf, we get the 
following linear relations between the components of the chiral matrix ps and 
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the on-shell functions fa, a = ±,1,2: 



Til 1 



777* 1 

Pk = 7^if- + n) + 7;ifl+^f2), (5.7) 

where ko = kQ — V{z) and km = (^q — "o^^)^^^- Formally these linear relations are 
invertible whenever fc^ 7^ 0. However, the on-shell momentum km may become 
imaginary, depending on the value of ko with respect to V. In this case, one needs 
to set /± = 0, since it actually corresponds to a situation where the mass-shell 
functions do not contribute in ^^-integration at all due to imaginary roots in the 
— m^)-factor. In general, the spectral approximation for the correlator gf, 
based on the mean field limit with vanishing gradients, is expected to become 
better when the distance from the potential wall increases. Asymptotically far 
away from the wall the spectral decomposition and thus the relations (15. 7p should 
be exact. Now, if we formally replace km ^ \km\ in Eq- (15.71) . we see that the 
resulting relations are invertible for all km ^ and there are no problems with 
imaginary km- Thus, we can always parametrize ps in terms of effective /- 
functions via these connection relations, which will reduce to the actual on-shell 
functions only asymptotically when km is real (or f± = 0). The use of such 
a parametrization is that the equations of motion for the effective /-functions, 
resulting from the matrix equation (15. 6p . become very simple inside the regions 
I and II (we set here mj = for simplicity): 

sdj^ = 
sdJl = ~2kof^ 

sdj^ = (5.8) 
ko 

Moreover, the structure of the dynamical matrix equation (15.61) readily implies 
that the only consistent condition for matching at 2; = is that all of its com- 
ponents are continuous, while the off-diagonals have kinks i.e. their derivatives 
have finite discontinuities. 



Real q, partial reflection 

Let us first consider the case with real q. The equations (15. 8p now imply that in- 
side the regions I and II the functions /|. are constants, so the asymptotic bound- 
ary conditions for the incoming fluxes fix f^^^ = 1 and = throughout the 
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regions. Moreover, the functions /i 2'^^^ are oscillatory in both regions. Hence, 
the asymptotic boundary condition for the coherence to vanish as z — >■ — cxd, 
kills it completely in region II: = 0- Using these conditions and the four 
matching conditions for the components of and (}} at z = 0, we can fix all the 
remaining constants: /+^\ as well as two constants related to the oscilla- 

tory coherence solutions in region I. We find that the reflected and transmitted 
fluxes are related by flux conservation: 

/f) = i^ and /im = l-/f' = (5.9) 

while the oscillatory coherence solution in region I is given by 

ff\z) = ^^cos{2kz) and ff\z) = -^sdjf\z) , (5.10) 

where we denote 

— ik — 2- (5-11) 

u!{UJ — V) — ni'' 

Imaginary q, total reflection 

In the case of imaginary q we have an additional boundary condition. That 
is, we cannot have any asymptotic mass-shell solution in region II: ff'^ = 
throughout (since they are constants). However, we see from Eq. (15.81) that the 
coherence solutions are exponentials in region II for imaginary q, so that the 
asymptotic boundary condition at z ^ —00 does not kill the exponentially de- 
caying mode. An otherwise similar analysis gives now the result with a complete 
reflection: 

/+^'^ = 1 and /5") = 0, (5.12) 
while the coherence solutions are 

with f2^^\z) = —sdzfi^^\z)/2uj in region I (there was a misprint in Eq. (118) in 
ref. [I]), and 

with f2^^^\z) = —sdzfi^^^\z)/2{uj — V) in region II. Note that the results in 
Eqs. (15.9115. 14l ). which agree with the standard Dirac equation approach [104], 
would not have been obtained, should we have dropped the coherence shell 
solutions /f 2 from our analysis. 



54 



5.1.2 Scalar bosons 



For scalar fields the free field KB-equation with a 2;-dependent potential becomes 

([fco - V{z)fe-^^^^^^'^ -kl- \dl + %Kd, - m')zA< = , (5.15) 

where we have taken = 0. Taking the real and imaginary parts of Eq. (I5.15p 
we find: 

(kl-\dl-e„,co^{\dXd^S)^^< = (5.16) 
(fcA + A;^sin(ia]'4j)2A< = 0, (5.17) 

where km{z) = ((/cq — V{z)y' — m?y/'^. We find again that the spectral phase 
space structure of the correlator A<, obtained by analyzing these equations in 
the mean field limit, is identical to the case of a spatially varying mass given by 
Eq. fl4.65p . except for the replacement kQ — V. We proceed now analogously 
to the spatially homogeneous case in section I4.4.2[ The n-th moments of the 
correlator A^ are defined as integrals over kz'- 

/(\k 
-^k:zA<{ko,k, = 0,kz,z). (5.18) 

By taking the 0th moment of Eq. fl5.16p and the 0th and 1st moments of 
Eq. fl5.17p we get the following closed set of equations for the three lowest mo- 
ments P0A,2- 

^dlpo + kl^po- P2 = 
dzpi = 

dzP2-\{dzkl)po = 0. (5.19) 

When the spectral correlator Eq. fl4.65p is used as an ansatz for A^, we get the 
following relations between the lowest moments po,i,2 and the on-shell functions 

f±,c- 

po = sgn(A;o-\/(^))-^(/+ + /-) + /. 

pi = sgn{ko-V{z))^-{U-f_) 

P2 = sgn{ko-V{z))^{U + U). (5.20) 

These linear relations are invertible whenever km 7^ 0, but for the case of imag- 
inary km we have to set /± = 0, as for fermions. Again, the spectral approxi- 
mation for the correlator, and consequently the relations f l5.20p . become exact 
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asymptotically far away from the wall. The matching conditions at 2; = in- 
duced by the moment equations (15.191) are now more complicated than with 
fermions. First, we note that dzPi vanishes everywhere, so pi is a constant 
throughout. The last Eq. fl5.19p implies that p2 must have a finite discontinu- 
ity over the barrier, which can be computed by integrating it over a step from 
z = —e to z = e to give p\ — = ^[k"^ — q^]poiz = 0). Using this in the first 
Eq. (15.191) . we see that also d'^po has at most a finite discontinuity over the step, 
implying finally that po and its derivative dzPo are continuous at z = 0. 



Real q, partial reflection 

We consider first the case with real q. The two latter Eqs. (I5.19P imply that 
dzPi,2 = 0, i.e. pi^2 are constants inside the regions I and II. These constants 
are partially fixed by the asymptotic boundary conditions /I = 1 and = 
through the relations (15.201) . Moreover, from the first Eq. (15.191) we now find 
that Po is oscillatory with a constant shift of amount P2/^m i^i both regions I 
and II, respectively. Combining Eqs. (I5.20p and (15.191) we then find that asymp- 
totically the coherence solutions are also oscillatory: fc = — 9^po/(4/c^), so the 
asymptotic vanishing of the coherence as ^ — — 00 kills the oscillatory part of 
Po completely in region II. Using these conditions and the four matching condi- 
tions for po,i,2 and dzPo at z = described above, we can fix all the remaining 
constants to obtain the solution 

I _ {k-qf rii_._ri_ 4A:g 

while the coherence solution in the region I is given by 

fl = l^/ficosi2kz). (5.22) 

Note that these solutions describe the shell structure of the correlator A*^ exactly 
only asymptotically as 2; ^ ±00. However, as in the case of fermions, replacing 
km ^ \km\ the relations are invertible for all km 7^ 0, so that the moments Po,i,2 
can be completely parametrized with these /-functions. In this (effective) sense 
the solutions (I5.2ip and (I5.22p are exact throughout the regions I and II. 



Imaginary q, total reflection 

As for fermions, in the case of imaginary q we cannot have any mass-shell so- 
lutions in region II, so that asymptotically both /j? = 0. Now the coherence 
solution fc = — (9fpo/(4A;^) is a superposition of decaying and increasing expo- 
nentials in region II, so that the asymptotic boundary condition at z —>■ —00 
does not kill the decaying mode. Because of these differences we get the solution 
with total reflection: 

/; = 1 and /i^ = 0, (5.23) 
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while the coherence solution is given by 



in region I, and 

= t4n^^'* (5-25) 

in region II. Again, we would not have obtained the results in Eqs. (15.21115. 25l) 
that agree with the standard Klein-Gordon approach, should we have dropped 
the coherence shell solutions fc from our analysis. 



5.2 Quantum reflection from a CP- violating 
mass wall 

Let us next consider the quantum reflection of free fermionic fields from a smooth 
CP- varying mass wall. This simple example is of relevance for electroweak 
baryogenesis (see refs. [105-107] and section [2. l.ip . where the complex spatially 
varying mass function arises from Yukawa couplings to the vacuum expectation 
value of the Higgs field(s): 

m{z) = y(j){z) , (5.26) 

where y is a Yukawa coupling and 0(^) is a complex scalar field corresponding 
to the total effect of the VEVs of (possibly) multiple Higgs fields [107]. For the 
complex mass wall we have used the following parametrization (see Fig. 15.21) : 

|0(z)| = i(l -tanh(z/C)) , arg[0(2)] = {l + tiin\i{z / 1^)) , (5.27) 

where £^ is the width of the wall, and is the total magnitude of the change 
of the phase of 0. The problem is described by the same integrated evolution 
equation as for the fermionic Klein problem, Eq. (15.61) . and the connection rela- 
tions Eq. (15. 7p between the components of the chiral matrix ps and the on-shell 
functions fa with the replacements V{z) —>■ and m^ij —>■ mRj{z). For the 
asymptotic boundary conditions we set the incoming flux from deep in the sym- 
metric phase at z — CX3 to unity and take no incoming flux from deep in the 
broken phase at z — > — oo, and correspondingly set the asymptotic coherence to 
zero as well: 

fl = l, z-^oo 
/; = /r = /2=0, z^-oo. (5.28) 

Let us summarize the results of the numerical calculations presented in ref. [I]. 
The calculations were performed using Eq. (15. 6p in Bloch representation: ps = 
^iido) + {9'^) ' together with the relations Eq. (15.70 and the boundary con- 
ditions Eq. (I5.28P . In figure 15.31 we plot the results for the case of s = 1 and 
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Figure 5.2: Reflection from a spatially varying mass wall, with the momenta of 
incoming (red) and outgoing (green) particles described by arrows. Spin s is 
conserved in the coUisionless case. 

q/\m^oo\ = 0.088, where q = (u;^ — Im-ooP)"*^^^ is the asymptotic momentum 
in the broken phase. In the left panel we show the values of the on-shell func- 
tions fa while in the right panel we plot the components of the chiral matrix 
p. We can see that in the symmetric phase (to the right from the wall) the 
system is a coherent superposition of left and right moving states with opposite 
fc^-momenta, and the /c^=0-shell functions are oscillating coherently. In the bro- 
ken phase however, all but the /--function die off and the state soon becomes 
a pure transmitted left moving state. Note that the physical flux is conserved: 
f+{oo) + /_(— oo) = /-(oo) = 1. For the chiral components we see that the 
imaginary part of plr goes to zero when z ^ — oo, as a result of our choice that 
m becomes asymptotically real in the broken phase. The diagonal components 
of p become large in the broken phase, because they represent chiral densities 
that are enhanced in the regime of small local velocity Vz = kz/oj due to flux 
conservation. 

In figure Em we show the particle-antiparticle flux-asymmetry Aj+ = (/+ — 
f^)z=oo as a function of g/|m_oo|. For antiparticles we simply need to make 
a replacement m — m* in all equations. The characteristic peaked shape of 
the flux-asymmetry can be understood as follows: The reflection amplitudes 
for both particles and antiparticles tend to unity in the limit of total reflection 
g 0, so the asymmetry must tend to zero. For large q on the other hand, the 
(anti)particles start to behave classically compared to the width (and height) 
of the wall, and hence the reflection amplitudes and the asymmetry start to 
decrease exponentially. 

These results have been derived earlier using the Dirac equation approach 
(see e.g. ref. [107]). Here we just wanted to demonstrate how the results are 
obtained using our formalism. We emphasize that for noninteracting fields our 
method is exact, as in the case of Klein problem. That is, all gradients of the 
mass are vanishing asymptotically at 2; ^ ±00, and consequently the mean 
field approximation becomes exact in that limit. So the /-t-functions indeed 
correspond to the physical fluxes asymptotically, and this is the only place where 
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Figure 5.3: a) Shown are the mass-shell functions f± corresponding to the left 
and right moving fluxes and functions /i_2, which encode the quantum coherence. 
We have taken s = 1 and g/|m_oo| = 0.088. b) The chiral density matrix 
components for the same solution. For the wall width and the total change of 
the phase we used = 2 and A9 = — 1. 




qAm(-oo)\ 

Figure 5.4: Shown is the current asymmetry of reflected states as a function 
of the asymptotic momentum to mass ratio in the broken phase. The wall 
parametrization is the same as in figure 15.31 
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this identification is needed in our analysis. Between the boundary regions the 
/-functions are just a parametrization of the chiral matrix p. However, the whole 
point of our approach is the ability to include the effects of collisions together 
with nonlocal coherence. With collisions the picture changes as we need the on- 
shell identifications of the /-functions for the reliable computation of the collision 
ternl^. For this reason our method becomes an approximation in the case with 
collisions, with the better quantitative results expected the thicker the wall is. 
But even in the thin wall limit we expect to get at least the correct qualitative 
behaviour, yet it would not be surprising if the results would prove to be quite 
accurate also quantitatively, since the region of (rapidly) varying mass is very 
narrow in this case, and outside the wall region the mean field approximation 
should work well. Moreover, the /-functions are presumably smooth also in the 
wall region, as suggested by the plots in Fig. 15.31 thus the overall error in the 
dynamics arising from the (approximate) computation of the collision term in 
that region should be rather small. 

5.3 Preheating of decaying fermions 

As a first example of our formalism including collisions, we consider the produc- 
tion of fermionic particles at preheating, during which the fermion is subjected 
to decays. A similar fermion production scenario in out-of-equilibrium condi- 
tions has been recently studied using the 2PI approach in ref. [108], where a 
complete numerical next-to- leading order calculation in a 1/A^-expansion of the 
2PI effective action was performed for a SU{2)l x 5'f/(2)ij ~ 0(4) symmetric 
theory with Yukawa interactions between fermions and bosons. As discussed 
in section 12.21 in preheating the (coupled) fermions obtain a rapidly oscillating 
effective mass due to the couplings to the inflaton condensate, for which we 
use here a simple cosine function (here t corresponds to a time coordinate in 
the conformal frame, and the expansion of the universe is neglected during the 
preheating) [75,102]: 



where tuq, A, B and uj^p (the inflaton oscillation frequency) are real constants. 
The decay of the fermion into daughter particles is modelled here by a left-chiral 
non-diagonal Yukawa interaction defined by the Lagrangian: 



where if) is the fermion field considered, q is some other fermion (quark) field 
and is a complex scalar field. For simplicity, we assume that the fields q and 
(f) form a thermal background i.e. they are in thermal equilibrium throughout 

^Also the corrections from mass gradients that arise from the gradient expansion of the 
coUision term become important 



m{t) = rriQ + A cos{2uJ^pt) + iB sin(2co'^t). 



(5.29) 



Ant = -y ^L(f) Qr + h.c. 



(5.30) 
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the preheating. A more reahstic but complex scenario could be modelled by 
taking into account the full dynamics of the daughter particles with appropriate 
equations of motion. By the assumption of thermal equilibrium for the daughter 
fields, the self energies T,^'^ for the fermion ip are related by the KMS relation 
E^(A;) = e^'^°S^(fc), and the interaction width F can be written in the form [II] 

r(A;) = i(l + e^'')tJ:<ik) = [ro7° - 7)] Pl , (5.31) 

where Pl = |(1 — 7^) is the left chiral projector and ro,3(A;) are real functions. 
Using these expressions and the mean field spectral function Eq. f l4.78p we find 
that the zeroth order gradient contribution to the collision term in the dy- 
namical equation (I4.66P is given by 

CH=-HD,9k-i9;:U} , (5.32) 

where 

D = ^{l + p,)Th with Th = To-hT^ (5.33) 

and (fi'^)eq is the thermal equilibrium limit of with fl^ = = l/(e^^'^fe + 1). 
Upon integration over /cq we get then the matrix equation (14.671) : 

dtPh = -t[H,Ph]-t{Ci^), (5.34) 
where the collision integral is now given by 



mO 



I fh _ fC<i\ _ l,k_( fh _ i-eqN J2L( fh _ f 

Uo Jo J "-(^Ua J3 J 2<^w3 J: 

7(/^-/r) 



m 
2uj 



- hV 



77l3 



( fh _ /eqx _ Uk ( fh _ feq\ m ( fh _ fee 

Us is J "-uiyjo Jo ) 2uj^Jo Jo 
C(/o'-/o') 



W^^^fi) j' ^'-''^ 

where we use the shorthand notations = f^^f^ and k = \k\. The Fj- 
functions appearing as coefficients are now projections on the different singular 
shells: 

rm(o,3)(|^M) = ro,3(fco = t^(t), 1^1) (positive mass-shell) 

Tood^l) = ro(A;o = 0,1^1) (A;o = -shell) . (5.36) 
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The explicit expressions for these F-functions are not presented here, but they 
can be read from the relations ro,3(A;) = |(1 + e^''°)iT,Q^{k) and the equations 
(6.23-6.28) in ref. [II]. Note that F on the negative mass-shell is related to that 
on the positive with ro,3(— fco, 1^1) = ^^o,3{ko, \k\), and furthermore that 
vanishes on the fco = -shell. 

Let us now summarize the results of the numerical calculations for particle 
production through a parametric resonance. A more detailed discussion is found 
in ref. [II]. In the calculations we have used the equations f l5.34p -( l5735l) with the 
connection relations fl4.69p . and the oscillatory mass function Eq. (15.291) . We 
have studied the time-evolution of the fermionic particle number and coherence 
with varying interaction strengths. In all cases, the initial condition for the 
evolution has been an uncorrected vacuum state with = ra^^ = /i^2 = 0. 
The results of our calculations for helicity h = —1 are presented in Fig. 15.51 In 
the upper panel we have considered the noninteracting case with F = 0. We 
see that the produced particle number (thick black line) as well as the "total 
amount" of coherence fc = \/ /i + /I (dotted blue line) increases steadily as a 
function of time. This increase takes place around the resonance peaks while 
between the peaks the particle number is essentially constant and the coherence 
oscillates with a constant amplitude, saturating to a maximum after a few res- 
onance crossings. This picture of parametric resonance depends sensitively on 
the parameters of the mass oscillation and the size of momentum \k\ and helic- 
ity h. For example for the opposite helicity h = +1 with otherwise the same 
parameters we would get a completely different figure without a clear resonance 
behavior. 

In the lower panel of Fig. 15.51 we have considered the case with interactions, 
F 7^ 0, with otherwise the same parametrization as in the noninteracting case. 
The difference compared with the noninteracting cases is quite dramatic. Now 
we see that the particle number drops between the resonance peaks, which results 
from the decays to the daughter particles, only to be regenerated again in the 
next resonance crossing. Also the growth of the coherence is now damped in 
comparison with the noninteracting case, and both the particle number and 
the coherence evolution settle into a stationary pattern after a few oscillation 
periods. In Fig. l5.6l we plot the particle number evolution for varying interaction 
strengths, again with otherwise the same parametrization as before. We see that 
the same pattern remains with increase in particle number during the resonance 
peaks, while the magnitude of damping (decays) between the peaks depends on 
the strength of the interactions, as expected. 

These results show that if the fermion has strong enough interactions with 
other fields during the preheating, the effects (especially decoherence) on the 
amount of produced particles can be significant. Caution must be exercised 
in that our approximation scheme with singular shell decompositions is not 
guaranteed to produce correct quantitative results in this example, because the 
rapidly oscillating mass function gives rise to big gradient corrections. However, 
at least the qualitative picture is expected to be correct. 



62 




Figure 5.5: Shown is the mean field number density n^^ of produced fermions, 
and the total amount of coherence fj! = {{fiY + (/2)^)^^^ (thin dotted blue 
line) for the negative helicity state h = —1. Effects of infiaton oscillations are 
modelled by a varying mass m(t) = (10 + 15 cos(2ci;^t)+i sin(2u;^t))T. The upper 
panel corresponds to the collisionless case, while in the lower panel the collision 



terms of Eq. fl5.35p were used, with the parameters \k\ = T, y = 5, 



o.o2r. 



= O.IT, and = T, where the temperature T sets the scale. Initially, at 
r = ujir,t = the fermion system is taken to be an uncorrelated vacuum state. 
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Figure 5.6: Shown is the number density in the same setting as before with 
changing interaction strengths. The thick black hne is the free field case with 
a coupling constant y = 0. The other lines are interacting cases, with y = 1 
(green line), y = 3 (blue dotted line) and y = 5 (red dashed line). 



5.4 Coherent production of decaying scalar 
particles 

In this section we consider a similar coherent particle production with decays 
as in the previous section 15.31 but now for scalar fields. A comprehensive study 
of resonant particle production using the 2PI approach has been carried out 
in ref. [109], where a complete numerical next-to- leading order calculation in a 
1/iV-expansion of the 2PI effective action was performed for an 0(A^)-symmetric 
scalar theory. We are not trying to model the scalar preheating scenario here, we 
simply want to compare the behaviour and the results of decaying scalar fields 
with that of fermions. For this purpose we take the mass term driving the scalar 
particle production to the absolute value of the fermionic mass in Eq. (15.291) : 

m^(t) = |mo + A cos{2uj^t) + iB sm{2uj^t)\'^ , (5.37) 

where mo. A, B and uj^ (oscillation frequency of the driving field (p) are real 
constants. For the decays we use the following Yukawa interaction: 

Ant = -?/^^0, (5.38) 

where is the real scalar field considered and ip is some fermion field, which is 
assumed to be in thermal equilibrium throughout the evolution. Hence the scalar 
field self energies II^''* are related by the KMS relation 11-^ (A;) = e^^°'n.'^{k) and 
the interaction width F is given by 

T{k) = ]^{l + e^^'')iYi<{k). (5.39) 
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Using these expressions and the mean field spectral function in Eq. fl4.80l) we find 
that up to first order gradients in the collision term iC^^n, the collision integrals 
in the dynamical equations fl4.73p are given by 

(koC-) = -^r„ [(/+-/.)- (/:^ - r)] , (5.40) 

where we have now defined = — 1), and we have neglected terms 

of ordeiEl OiV"^ ^Tdtrm?). The Fj-functions appearing in Eq. (15.401) are again 
projections onto the mass- and the coherence shells: 

Trnd^M) = '^{ko = oj^{t),\k\) (positive mass-shell) 

9fc„ro(|fc|) = 9fc„r(fco = 0,|fc|) (fco = -shell) . (5.41) 

As in the previous section, we present no explicit expressions for these Fj- 
functions here, but they can be read from the equations (15.391) and (A.5-A.11) 
in ref. [III]. Note that F on the negative mass-shell is again simply related to 
that on the positive mass-shell: F(— fco, 1^1) = — F(/co, \k\) and it vanishes on 
fco = 0-shell. This is the reason that we need to include the first order gradients 
of the collision term in the scalar case, as we see that dk^Vf) is the lowest order 
nonvanishing fco = 0-shell contribution, giving rise to decoherence effects. Using 
the relations (14.741) we can now express the on-shell functions f±^c in terms of 
the moments Po,i,2 to find a closed set of equations for these moment functions: 

dtPl = —Tm (Pl - Pl,eq) 

^k 

dtp2 - ^dt{m^)pQ = F„ (p2 - P2,eq) , (5.42) 

^ ^k 

where pj eq are the the thermal equilibrium values for the moments with f± = 
and fc = 0. 

Let us now summarize the results of the numerical calculations for particle 
production through a parametric resonance. A more detailed discussion can 
be found in ref. [II]. It turns out that the equations of motion (15.421) become 
numerically unstable for rapidly oscillating driving mass terms in Eq. (15.371) . For 
that reason we have written Eq. (I5.42p in a different form with nonlinear terms 
and an integration constant, given by Eqs. (8.3)-(8.5) in ref. [Ill], for which 
the stability problems do not occur. Using these equations we have performed 



^ These terms are neglected also indirectly, by recursive use of equations of motion (|4.73p 
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Figure 5.7: Shown is the number density (thick sohd hne) and the coherence 
function fc{\k\) (dotted(blue) hne) with changing interactions. The driving mass 
function is taken to be m^(t) = |(1 + 1.5 cos(2c<j^t) + i 0.1 sin(2ct;^t))Tp. The 
upper panel corresponds to the case without colhsions. In the central panel we 
have included the collision terms on the mass-shells, but kept dk^To = 0. In the 
lowest panel the full interaction terms were kept for all shells. For parameters we 
have used \k\ = 0.6 T, y = 1, = 0.1 T and = O.IT, where the temperature 
T sets the scale. Initially, at r = uj^t = the system is in the adiabatic vacuum. 
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Figure 5.8: The same configuration as in the lowest panel of Fig. 15.71 but with 
the driving mass term smoothly set to a constant m = 2.5 T after r > Att. 



similar calculations as for fermions in section [531 z.e. we have studied the time- 
evolution of the particle number and coherence from an uncorrelated vacuum 
state rij^ = Uj^ = fc = 0, with varying interaction strengths. 

The results of the calculations are presented in Fig. 15.71 In the upper panel we 
have considered the noninteracting case with F = 0. The increase of the number 
density (thick solid line) is seen to be accompanied by a steady growth of the 
amplitude of the coherence (thin dotted blue line). Similar to fermions there 
are regions (resonance "peaks") where the increase takes place, while between 
the peaks the particle number and the amplitude of the coherence are relatively 
constant. In the middle panel we considered the case where only the mass- 
shell collision terms were included but we set artificially dko^Q = 0, while in 
the lowest panel all collision terms were included properly. In these cases with 
interactions we see again that the decrease of the particle number (decays to 
daughter particles) and of the amplitude of the coherence (decoherence) takes 
place between the resonance peaks. By comparing the two lowest panels we can 
conclude that the dominant effect of interactions comes from the coherence shell 
contribution, dk^To, which would be absent in the standard quantum Boltzmann 
approach. 

To see the effects of thermalization and decoherence more clearly, we consider 
finally a case where the driving term was switched off smoothly at a specific mo- 
ment. This scenario is depicted in Fig. 15.81 We can see a smooth thermalization 
of the particle number towards the equilibrium value, accompanied by decoher- 
ence towards zero in the form of damped oscillations of the coherence solution 
fc- From the master equations fl5.42p we can identify that the corresponding 
thermalization and decoherence rates are given by F^ and dk^To, respectively. 
This example of thermalization and decoherence of a highly coherent initial state 
provides a more suited application for our eQPA scheme, as the big gradients of 
the background field are no longer a problem. It would be interesting to study a 
more realistic example of scalar field thermalization using our methods, e.g. with 
a model of quartic self interaction. 
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5.4.1 Resummation of the coherence contributions in 
the coUision term 



In the examples of sections 15.31 and 15.41 the gradient expansion of the coUision 
term was truncated to the lowest order that includes decoherence. For fermions 
this was the zeroth order and for scalars the first order. However, as discussed 
in section IT7t also the higher order coherence contributions in the collision term 
should be taken into account, since their expansion is not (completely) controlled 
by small parameters. In this section we present an analysis that shows how the 
coherence contributions in the collision term can be resummed for a scalar field 
with a thermal interaction. The results of this section were not presented in 
refs. [II, III], and a new publication [110] on these and the corresponding results 
for fermions is being prepared. 

To be specific, we want to expand the collision term to leading order in F 
(and dtm?): Oi = 0{T, dtm^), neglecting the terms of order O2 = 0(F^, F dtm?) 
and higher. We begin by an observation that the on-shell functions obey the 
zeroth order equations: 

dtf± = Ci (5.43) 
d^J^ + Aulf, = Oi, (5.44) 

which follow directly from Eqs. fl4.73p and the relations fl4.74p . Furthermore, 
the first part of the (standard) collision term (13.441) gives 

e-**{F}{zA<} = exp )f2A< + C2 , (5.45) 

since dtV is proportional to dtm^ (as well as to F) for the interaction with thermal 
background. When the fco-integration is performed for a term of order n > 1 of 
this expansion, we see that the only leading order contribution is the projection 
to the fco-shell: 

/ ^ ^r^^< = dmk.=o dUc + O2 , (5.46) 

where we have used Eqs. (15. 43 p . Now, using Eq. (15.44p we find that to zeroth 
order all the derivatives d^fc are proportional to either or dtfc- d'^'^fc = 
{-AulYfc + Oi and S^^+V, = {-AulYdtfc + Ci- Using these relations we get 
then 

e~**{F}{iA<} = cosh(a;^afejF|fco=o/c + «— sinh(a;^afcjF|fco=o^i/c 
+regular + O2 , (5.47) 

where regular refers to standard first order mass-shell contribution. The second 
part of the collision term (13.440 . proportional to the spectral function A, gives 
only the regular contribution, since the spectral function does not involve the 
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coherence contributions in our approach. Putting things together, we find the 
following generalizations of Eq. (15.421) : 



^d'^po + uj^po - p2 



2uo, 



sinh(u;g9fcjr|fco=oc^tPo 



dtPi = — —Tm (pi - Pi,cq) - cosh.{uj^dko)T\ko=o(po - ^^2) 



dtp2 - -dt{rn?)pQ 



-r.m (P2 - P2. 



(5.48) 



where the terms of order O2 are neglected. We see that the infinite /co-derivative 
expansions of F in Eq. (15.481) are not controlled by the small parameters F and 
dtrr?\ instead they provide expansions in powers of (ioj^ or cj^/l/cl etc. that 
are typically of order unity. Moreover, because of the generally nontrivial /cq- 
dependence of the (thermal) F, a brute force calculation of the (infinitely many) 
derivatives seems impossible. However, it turns out that these expansions can 
be resummed by a neat trick. That is, by denoting F(ro) the Fourier transform 
of F with respect to /cq and using the identity sinh(ia;) = 'isin(x), we find 



sinh(a;^9fco)F|fco=o 



dro sinh(iu;^ro)F(ro)e' 
1 



fco=0 



dro^(e*-s^° 



)r(ro) 



and similarly 



cosh(a;^(9fco)F 



1 



F(-w. 



fco=0 



[F(.;,^) + F(-^,^)]. 



(5.49) 
(5.50) 

. -r(^,^), 

so that we get simply: sinh(u;^9fco)F|fco=o = r(u;^) = F^ and cosh(a;^9fcjF|A:o=o = 
0. When these amazingly simple results are substituted in Eq. (I5.48P we get the 
final equations to the leading order as 



For the thermal interaction considered we furthermore have: F( 



-a^Vo + ^\pQ - p2 

dtpi 



dtP2 - ^dt{m^)pQ 



2uo, 



-Fm dtPQ 



-Fm (pi - Pl,eq) 



-Fm (p2 - P2, 



eqy ; 



(5.51) 



where the terms of order 0{T^ ,T dtm^) and higher have been consistently ne- 
glected. We see that due to the resummation of the coherence contributions, the 
resulting equations have actually simplified, as F^ is now the only interaction 
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T 

Figure 5.9: The number density (thick sohd hne) and the coherence function 
fc{\k\) (dotted(blue) hne) in the cases of non-resummed (upper panel) and re- 
summed (lower panel) coherence shell collision terms. The configuration is the 
same as in the lowest panel of Fig. 15.71 

term to be evaluated. However, what is more important, these equations provide 
a controlled expansion of the collision term to the leading order in F and dtw?. 

Let us now compare the numerical results obtained by using the resummed 
equations f l5.5ip to the original ones, presented in Figs. 15. 7115. 81 In figure [53] we 
show the results with the original (non-resummed) collision term in the upper 
panel and the resummed one in the lower panel. We have used the same con- 
figuration as in Fig. 15. 7[ so the upper panel here is just a copy of the lowest 
panel in Fig. 15. 7[ We see only a slight difference in the results, with the net 
amount of produced particles slightly increased in the resummed case. This is 
explained in that the average value of the oscillating resummed decoherence rate 
Tm/oj is approximately equal to the constant non-resummed rate dk^TQ in this 
particular case. In the example of thermalization and decoherence, however, a 
bigger difference can be seen. The case with the resummed collision term and 
otherwise the same configuration as in Fig. 15.81 is presented in Fig. 15.101 We 
see that decoherence effects are clearly stronger in the resummed case, because 
now Tm/oj freezes to a value that is significantly larger than SfcgFo. The ther- 
malization rate of the particle number, on the other hand, is not changed at all. 



70 




0.11 

% 8 10 12 14 16 18 20 



Figure 5.10: The same configuration as in Fig. 15.81 but with a resummed coher- 
ence shell collision term. 

because the particle number and the coherence solution decouple in the constant 
mass limit in the case of a simple thermal interaction, and thus the (increased) 
decoherence rate does not affect the thermalization of the particle number. 
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Chapter 6 

Conclusions and outlook 



We have developed a novel approximation scheme (eQPA) that enables us to 
include the effects of nonlocal quantum coherence in the standard kinetic ap- 
proach to nonequilibrium quantum dynamics. The key element in our scheme is 
the finding of new singular shell solutions in the phase space of 2-point Wight- 
man function, located at = for spatially homogeneous problems and at 
kz = for a static planar symmetric case. We have interpreted these new so- 
lutions as describing the nonlocal coherence between the "opposite" mass-shell 
excitations with either ko = ±a;^ or k^ = ^km, respectively. This nontrivial 
phase space structure is then inserted in the dynamical Kadanoff-Baym equa- 
tion for the 2-point correlator, leading to a closed set of transport equations for 
the corresponding on-shell distribution functions /, that provides an extension 
of the standard quantum Boltzmann equation to include nonlocal coherence. 

We have shown how the eQPA scheme emerges from first principles in chap- 
ters [HMl We used the 2PI effective action method to derive the self-consistent 
Schwinger-Dyson equations for the 2-point correlation functions of the system. 
We then performed a partial Fourier transformation to obtain the Kadanoff- 
Baym transport equations for the Wightman function in the mixed representa- 
tion. When these equations were analyzed in detail in the kinetic regime with 
weak interactions and a slowly varying background field, the nontrivial singular 
shell picture with the new coherence solutions was found. 

The analysis for fermions was originally presented in refs. [I, II], and for 
scalar fields in ref . [Ill] . In these papers we have considered several applications 
to demonstrate the use of our formalism. The (coUisionless) Klein problem was 
considered in refs. [I, III] for fermions and scalars, respectively. In ref. [I] we also 
considered the (fermionic) quantum reflection from a CP-violating mass wall, 
and a topic omitted in the introductory part of this thesis, a generalization of 
our formalism to flavour mixing. The coherent production of decaying fermionic 
and scalar particles was considered in refs. [II, III] , respectively. Finally, the 
nonrelativistic limit of our formalism was considered in ref. [III]. 

In the coUisionless examples of the Klein problem and the fermionic quan- 
tum reflection from a smooth mass wall, presented in sections 15.11 and 15.21 our 
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formalism was seen to reproduce the same (exact) results as the standard ap- 
proach using the Dirac and Klein-Gordon equations. Of course, the whole point 
of our approach is the ability to include the effects of collisions together with the 
nonlocal coherence. We studied the role of collisions in the examples of coherent 
production of decaying particles in sections 15.31 and We saw a clear decrease 
in the amount of produced particles due to smooth thermalization and deco- 
herence effects caused by the collisions. A caution must be given, in that the 
parameters of the models were chosen such that the gradients of the background 
field were not necessarily small, and thus the numerical results of these partic- 
ular examples are perhaps only qualitatively correct. In section 15.4.11 we have 
shown how a resummation can be carried out for the coherence contribution of 
the (infinite) gradient expansion of the collision term. This procedure leads to a 
controlled leading order result in powers of the interaction width F and gradients 
of the background field. These new results are not presented in refs. [I, II, III] 
but will be published elsewhere [110], along with the corresponding results for 
fermions. 

One of the most interesting future applications for our eQPA scheme is a 
realistic calculation of the baryon asymmetry in electroweak baryogenesis, where 
the nonequilibrium dynamics of fermions in the neighbourhood of an expanding 
CP-violating phase transition wall needs to be considered. As discussed above, 
the relevant quantum reflection effects from the wall have not yet been considered 
in a consistent framework based on quantum field theory, including the effects of 
decohering collisions. Our eQPA scheme provides a tool for the study of such a 
problem with a smoothly varying background field and well-defined asymptotic 
regions. Indeed, the = 0-shell coherence solutions are just the missing piece 
for describing the nonlocal quantum reflection effects. In the case of a thick wall 
with small gradients, the mean field approximation works well throughout, and 
in the opposite limit of a thin wall at least the correct qualitative behaviour 
is expected. However, the results in this case should be quite accurate also 
quantitatively, since the region of (rapidly) varying mass is very narrow, and 
outside the wall region the mean field approximation works well again. 

Other interesting applications for our formalism include neutrino flavour os- 
cillations in inhomogeneous background, where the relevant coherence effects 
occur between different flavour states travelling in the same direction. More 
generally, we expect that a number of other problems that can be studied in 
the standard kinetic approach with the Boltzmann equations, could be studied 
using our eQPA scheme including the coherence effects. 
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